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We compute the dynamics of locahzed excitations produced by a quantum quench in the spin 
1/2 XXZ chain. Using numerics combining the density matrix renormalization group and exact 
time evolution, as well as analytical arguments, we show that fractionalization due to interactions 
in the pre-quench state gives rise to "ultrarelativistic" density waves that travel at the maximum 
band velocity. The system is initially prepared in the ground state of the chain within the gapless 
XY phase, which admits a Luttinger liquid (LL) description at low energies and long wavelengths. 
The Hamiltonian is then suddenly quenched to a band insulator, after which the chain evolves 
unitarily. Through the gapped dispersion of the insulator spectrum, the post-quench dynamics 
serve as a "velocity microscope," revealing initial state particle correlations via space time density 
propagation. We show that the ultrarelativistic wave production is tied to the particular way in 
which fractionalization evades Pauli-blocking in the zero-temperature initial LL state. 

PACS numbers: 71.10.Pm, 05.45. Yv, 64.60.Ht, 67.85.-d 



I. INTRODUCTION 

In the labyrinth of one-dimcnsional (ID) quantum 
many-body physics, the Luttinger liquid (LL) lurks 
around nearly every corner. It emerges as the low energy 
field theory description of interacting Bose gases, gapless 
quantum spin chains, fermion lattice models (Hubbard, 
etc.), electrons in metallic carbon nanotubes, and chiral 
quantum Hall edge statesii"— Luttinger liquid physics is 
universal: it reduces the complexities of myriad micro- 
scopic models to the hydrodynamics of free bosons. 

Despite its apparent simplicity, the LL description of 
interacting fermions exhibits a number of rather peculiar 
properties, due to the advent of quasiparticle fractional- 
ization. The elementary excitations of a LL are collective 
density waves that carry fractional (electric or number) 
charge, relative to the "bare" fermionic constituents; in- 
jecting a bare fermion into a LL causes it to "break up" 
into many pieces. This collectivization of the dynam- 
ics due to fractionalization leads to a host of predicted 
anomalies, including the low bias suppression of the tun- 
neling density of states, and perfect insulating behavior 
at zero temperature due to the presence of even a single 
impurityi^ For spinful fermions, fractionalization induces 
spin-charge separation.i^ Interestingly enough, zero tem- 
perature dc transport in a clean quantum wire through 
ideal Fermi liquid leads shows no signature of fractional- 
ization; the conductance is quantized to e^/ft- per channel, 
irrespective of the interactionsi^ 

In this paper, we describe a "transport" effect that 
directly exhibits fractionalization in a LL, observed in 
the dynamics of a density fluctuation following a sud- 
den quantum quench. In a quantum quench, a system 
is prepared in an eigenstate of some initial Hamiltonian. 
In our case, we take the ground state of an XXZ chain 
with a non-uniform density profile, which possesses a low- 



energy LL description. At the time of the quench, by 
external means a sudden deformation is affected upon 
the Hamiltonian, which subsequently drives the unitary 
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FIG. 1. Ultrarelativistic wave generated from an initial den- 
sity bump, following an interacting quench. The Luttinger liq- 
uid ground state of an interacting XXZ chain is time-evolved 
according to a non-interacting, band insulator Hamiltonian. 
The density 5p{t,Xj) due to the inhomogeneity is plotted at 
time slices t = 0, 12, 24, 36, and 48 after the quench; fainter 
(bolder) traces depict earlier (later) times. The evolution is 
symmetric about Xj — 0. In this figure, red dashed lines 
were obtained from a combination of DMRG and exact time 
evolution for the XXZ chain, while blue solid lines are the 
prediction of continuum sine-Gordon field theory. The curves 
marked "asymptotic" are the analytical result for the "reg- 
ularized supersoliton," obtained in Eq. p.48|) . The initial 
coupling strength is 7 = —0.872, corresponding to a = 0.7. 
The initial bump has width A = 12 and weight Q — 0.10; 
the mass gap is i\f = 1/8. The continuum data obtains from 
numerical integration of Eq. (|3.47p with a = 0.75 and C = 1- 
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post-quench dynamics. Here, the post-quench spectrum 
consists of non-interacting fermions, with a Hamiltonian 
that possesses a gapped, band insulator ground state (the 
XX chain in the presence of a sublattice staggered ex- 
ternal field). We show that the relative fractionaliza- 
tion of the pre-quench system (due to interparticle in- 
teractions) leads to the production of "ultrarelativistic" 
density waves after the quench. These waves travel at 
the maximum band velocity, and exhibit a particular 
shape set by the interaction strength. The propagating 
density waves are "elementary excitations" of the post- 
quench non-equilibrium state; they occur because the 
fractionalized density inhomogeneity of the initial LL "in- 
jects" high momentum excitations into the post-quench 
band insulator. By contrast, under the same conditions 
a quench from the ground state of the non-interacting 
Fermi gas (the XX point of the XXZ chain) into the band 
insulator yields only dispersive density dynamics, a con- 
sequence of Pauli-blocking. Our setup can be viewed as a 
"quench spectroscopy" of fractionalization in a Luttinger 
liquid. 
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FIG. 2. Dispersive decay of a density bump following a non- 
interacting (7 = a = 0), non-relativistic {MA = 3/2) quench. 
This is the same as Fig. [1] but for a quench from a non- 
interacting Fermi gas to a band insulator. Parameters A, Q, 
and M are as in Fig. [T] Red dashed lines are the results of 
exact diagonalization of the lattice Hamiltonian, while blue 
solid lines are the continuum predictions. 



In the last decade, rapid experimental progress^ in ul- 
tracold atoms and optical lattice gases has transported 
far-from-equilibrium many body physics fully into the 
quantum realm. In these systems, the quantum quench 
has emerged as a primary tool with which to investigate 
dynamics. Quenches have been performed in boson^i^— 
and fermioni^"— systems, with and without optical lat- 
tices, in one, two, and three dimensions. An ultra- 
cold gas can be very well isolated from its environ- 
ment, and provides an unprecedented degree of control 
in terms of realizing model systems and manipulating 
their parametersii Theoretical work has focused primar- 
ily on thermalization^^— , quantum critical scaling and 
defect production^"— and correlation functions in spa- 
tially homogeneous systems.—"— Prior art on Luttinger 
liquid, sine-Gordon, and XXZ chain quenches includes 
that of Refs. [13, [11,111 (H-IH MM- Wavepackets have 
been previously employed in the study of excitations in- 
duced by a local quench)^"— in which the Hamiltonian 
deformation is restricted to a spatial subregion of the 
larger system. The characterization of spatially inhomo- 
geneous dynamics following a global parameter quench 
(as studied here) is a more recent development)^"— i^i^ 



Many of the previous schemes proposed or executed in 
the theoreticali^i^iSad^d^ and experimental^i^iHiii liter- 
ature can be termed "hard quenches." In these works, 
large changes in a parameter value or trap geometry lead 
to the excitation of novel high energy states^'^ i"'^^'^^ whose 
physics has little to do with the low-energy sectors of ei- 
ther the initial or final Hamiltonians (for an interesting 
exception, see Ref . l47h . Our goal in this paper is different: 
we use a "soft quench" (defined below) as a low-energy 
probe of the initial state. 



A. Overview 

1. XXZ quench protocol; velocity microscope 

We study the dynamics following a quantum quench 
in the ID spin 1/2 XXZ chain. Working in the equiva- 
lent spinless (or spin-polarized) fermion representation, 
we investigate the time evolution of the particle den- 
sity induced by a non-uniform initial state. Other works 
treating XXZ and sine-Gordon quenches subject to ini- 
tial state inhomogeneity include Refs. l37l - [39t see Sec. IV Cl 
for a discussion. 

We consider a system initially prepared in the ground 
state of the XXZ chain in its gapless XY phase, subject to 
an external field. The field induces a localized "bump" in 
the density profile of the otherwise spatially homogeneous 
system. This state is further characterized by the spin 
anisotropy 7 of the pre-quench Sf_^_i coupling, i.e. the 
four fermion interaction strength. The gapless XY phase 
of the XXZ chain admits a low energy Luttinger liquid 
(LL) description.—!^ At time t = 0, the system Hamilto- 
nian is deformed discontinuously: 7 is set to zero, while a 
sublatticc-staggered external field is simultaneously ap- 
plied along the length of the chain, opening up a gap 
in the spectrum. In the fermion language, the ground 
state of the post-quench ( "final" ) Hamiltonian is a non- 
interacting band insulator with a doubled unit cell. 

The lattice quench with 7 = in the initial XY state is 
special, because both the initial and final Hamiltonians 
are non-interacting in the fermion language. We dub this 
the "non-interacting" quench; the exact solution can be 
written for the time evolution of the density expectation 
value. By contrast, for 7 ^ ("interacting" quench) the 
initial Hamiltonian is interacting in the fermion language 
and not soluble by elementary means. Although the XXZ 
chain is intcgrable, the non-uniform density profile makes 
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difficult the application of the Bethe ansatz method. In- 
stead, in this paper we use the Density Matrix Renormal- 
ization Group (DMRG) to numerically compute correla- 
tion functions of the initial ground state. For both the 
non-interacting and interacting quenches, the dynamics 
generated by the non-interacting band insulator Hamilto- 
nian are determined exactly. This allows us to avoid the 
use of more computationally intensive, time-dependent 
DMRG calculations. We exploit this advantage to ana- 
lyze larger system sizes than previous numerical quench 
studies of the XXZ chain^iM 

The idea behind this setup is to use the quench into 
a gapped, dispersive phase as a "velocity microscope" 
on the initial correlated LL state. The non-uniform ini- 
tial density profile creates additional excitations on top 
of the homogeneous bath induced by the global param- 
eter quench, leading to real space dynamics that can in 
principle be directly observed. Particles composing space 
time density fluctuations are excited with a broad range 
of momenta; these are velocity-resolved by the dispersive 
post-quench spectrum. By contrast, time evolution with 
a generic gapless post-quench Hamiltonian in ID (such as 
that governing a continuum conformal field theory) pro- 
duces only pure left- and right-moving "ultrarclativistic" 
waves, regardless of the structure of the initial state i^^i^^ 

Throughout this work we make the crucial assumption 
of a "soft quench," defined as follows. The magnitude 
of the gap in the post-quench Hamiltonian is specified 
by a dimensionless parameter Ma, where 1/Af gives the 
"Compton wavelength" for the low energy, massive exci- 
tations of the band insulator, and a denotes the lattice 
spacing. In addition, we assume a Gaussian density in- 
homogeneity in the initial state of width A. The assump- 
tion of a "soft" quench requires that 



««^<A, 



(1.1) 



i.e. that the low-energy Compton wavelength dwarfs the 
lattice spacing, while the width of the initial state inho- 
mogeneity exceeds the Compton wavelength. The first 
assumption guarantees that the gap opens in the low- 
energy sector of the band Hamiltonian. The second 
MA > 1 assures that any excitation of large-momentum 
particles post-quench arises from the correlated charac- 
ter of the LL, and not the excessive "squeezing" of the 
initial density bump. 

Despite the requirement in Eq. (|l.ip . we will consider 
quenches with "intermediate" to "large" values of the 
initial XXZ interaction strength 7, approaching the fer- 
romagnetic transition at 7 = — 1. It is far from obvious 
that a change from I7I ^ 1 to 7 = preserves the no- 
tion of a "soft quench" as articulated above. Because the 
low-energy description throughout the gapless phase is a 
LL, it is nevertheless the case. 
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FIG. 3. Post-quench evolution of a density bump: depen- 
dence on initial state interparticle interaction strength. Each 
subpanel exhibits a three-dimensional view of a lattice quench 
(obtained by DMRG -I- exact time evolution) into the gapped 
band insulator; Q, A, and M are the same as in Figs. [1] and [2] 
The four frames depict quenches with increasing initial state 
interactions, cr = (non-interacting), 0.4, 0.7, and 1.0 (top to 
bottom). The cyan line demarks the maximum band propa- 
gation velocity ("speed of light"), Umax(M = 1/8) ~ 1.77. 
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2. Sine-Gordon and "Supersolitons" 

The XXZ quench can be interpreted as a lattice version 
of the continuum sine-Gordon field theory analyzed pre- 
viously in Ref. l43l In that work, spatiotemporal dynam- 
ics were computed in a quench across a quantum critical 
point. In Ref. a LL ground state subject to an inho- 
mogeneous density modulation was time-evolved accord- 
ing to a translationally invariant, post-quench Hamilto- 
nian favoring a gapped Mott ground state. The Mott 
Hamiltonian generating the dynamics was chosen to re- 
side at the Luther-Emery^ point of the sine- Gordon 
model, where the excitation spectrum consists entirely 
of non-interacting, massive Dirac fermionsJ^i^ In the 
XXZ chain quench studied here, the final state band in- 
sulator carriers play the role of the non-interacting Dirac 
fcrmions that compose the sine-Gordon spectrum at the 
Luther-Emery point; the Mott gap of the sine-Gordon 
model is here substituted by the band gap. 

A localized density inhomogeneity in the sine- Gordon 
quench launches ultrarelativistic, non-dispersing travel- 
ing waves, dubbed "supersolitons" in Ref. |43. The super- 
soliton exhibits a rigid shape, propagates at the "speed 
of light" vf (the Fermi velocity), and possesses an am- 
plitude that grows in time as i"'/^. The exponent a > 
characterizes the fractionalization of the initial LL state 
relative to the final Mott insulator. For the case a = 
(non- interacting quench), there is no fractionalization 
and no supcrsoliton; the density dynamics of such a 
quench with MA > 1 show only dispersive broadening. 

In this work, we demonstrate that the supersoliton 
arises in the sine-Gordon quench for a > due to the par- 
ticular way in which LL fractionalization evades Pauli- 
blocking. This is made explicit through a calculation of 
the local phase space (Wigner) distribution in the pre- 
quench LL. The result is a power-law occupation of mo- 
mentum states in the post-quench insulator that trans- 
lates into a singular peak at vp the corresponding (lo- 
cal) velocity distribution. Because velocity is conserved 
by the post-quench Hamiltonian, the spectral weight as- 
sociated to the singularity is translated at vp. By con- 
trast, a non-interacting quench with cr = and MA ^ 1 
excites only small velocities v < vp/MA. 

A key point is that it is the long-distance behavior of 
correlations in the initial state that permits the evasion 
of Pauli-blocking in the fractionalized case. Although 
lattice details can and do modify the ultraviolet behav- 
ior of correlations in the XXZ chain considered here, the 
fundamental distinction between non-interacting and in- 
teracting quenches remains a robust feature of the soft 
quench satisfying Eq. (|l.ip . 

3. Preview of numerical results 

We defer a detailed discussion of our XXZ chain quench 
results to the main text; the impatient reader may 
consult Sec. |V] for a summary. Instead, we exhibit a 



few graphs that demonstrate the qualitative difference 
between the interacting and non-interacting quenches. 
Fig.[T]shows the "ultrarelativistic" density wave launched 
in an interacting quench satisfying the constraint in 
Eq. (|l.ip . (We set the lattice spacing a = 1). By con- 
trast. Fig. [5] depicts a non-interacting quench; in this 
case, only dispersive broadening of the initial density in- 
homogeneity is seen. The parameters in these two figures 
are the same, except for the interaction strength, quan- 
tified by a parameter 17(7). For a non-interacting quench 
one has 17(7 = 0) = 0; otherwise cr > and increases 
monotonically with |7|. The evolution of an XXZ chain 
quench as a function of the interaction strength a is de- 
picted as a 3D plot sequence in Fig. |31 

The blue continuous curves in Figs. [T] and [5] are ob- 
tained using an ultraviolet-regularized version of the sine- 
Gordon quench studied in Ref. |43| . The regularization 
models the effects of neglected lattice scale details in a 
very crude way. The regularization cuts off the amplifi- 
cation of the supersoliton predicted for the pure sine- 
Gordon quench; it also leads to a modification of its 
interaction-dependent shape. Because of the close agree- 
ment between the field theory and lattice results, we in- 
terpret the ultrarelativistic density wave appearing in the 
interacting quench (Figs. [T] and [3)) as a "regularized" su- 
persoliton. 

We emphasize that the quench dynamics described in 
this paper are fully quantum coherent; the absence of in- 
terparticle scattering in the post-quench band insulator 
prevents dephasing or thermalization. The "fractional- 
ized" density dynamics reflect the many-body entangle- 
ment of the initial gapless state. Future work incorpo- 
rating integrability-preserving interactions post-quench 
could prove particularly interesting, as discussed in the 
Conclusion. 



B. Outline 

The organization of this paper is as follows. In Sec. |lll 
we define the pre- and post-quench XXZ Hamiltonians 
and set up the dynamics to be computed. In Sec. IIIIl 
we provide a comprehensive analysis linking the XXZ 
chain quench studied here to the corresponding version 
in the continuum, low-energy sine-Gordon field theory. 
We begin in Sec. IIII Al with a quick review of single 
particle relativistic wavepacket mechanics, where we em- 
phasize the distinction between "relativistic" and "non- 
rclativistic" wavepacket propagation. In Sec. IIII Bl we 
describe the solution to the pure sine-Gordon quench. We 
identify the supersoliton, discussed previously in Ref. IH. 
The global and local (Wigner) distributions induced in 
the lattice and continuum quenches are discussed in 
Sec. IIII C[ wherein the origin of the supersoliton is re- 
vealed. In Sec. IIIIDI the ultraviolet modifications of the 
sine-Gordon theory necessary to model the lattice quench 
are articulated, and relevant time scales are defined. 

Numerical results obtained for the time evolution of 
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the XXZ chain quench are presented and discussed in 
Sec. HYl Results for the non-interacting and interact- 
ing quenches arc exhibited and compared to the regu- 
larized sine-Gordon theory. We summarize our conclu- 
sions in Sec. |Vl and finish with a discussion of open 
questions. The asymptotic analysis method used to ob- 
tain key analytical results is explicated in Appendix [X] 
Appendix |B] recapitulates the notion of fractionaliza- 
tion in a Luttinger liquid. In Appendix \C\ we derive 
the local (Wigner) velocity distributions induced by the 
initial state inhomogeneity, in the interacting and non- 
interacting continuum quenches. 



II. QUENCH SETUP 

A. Lattice model 

In a (sudden) quantum quench, one prepares the sys- 
tem in an eigenstate of an initial Hamiltonian , and 
subsequently time evolves under a different final Hamil- 
tonian, H^f\ We consider the XXZ spin 1/2 Heisenberg 
chain. 



(2.1) 



Via the Jordan- Wigner transformation, the spin chain 
is equivalent to a model of spinless (or spin-polarized) 
fermions whose Hamiltonian is given by 



H.c. 



(2.2) 

where J denotes the nearest-neighbor hopping ampli- 
tude, 2 J7 is a nearest-neighbor density-density interac- 
tion strength, and Hi represents a site-dependent chemi- 
cal potential. In Eq. (|2.2p . Ci and satisfy c,;cj + cljCi = 



5ij, and Sn 



cUi 



1/2. We will quench from a 
ground state in the gapless XY phase of this Hamilto- 
nian (labeled by the dimensionless interaction strength 
7) to a non-interacting, band insulator state. The latter 
is induced via the application of a unit cell doubling, 
sublattice-staggered chemical potential. Furthermore, 
by applying a localized Gaussian chemical potential, we 
will induce a density inhomogeneity into the initial state 
whose dynamics will reveal the effects of our quantum 
quench. 

The zero temperature phase diagram for the XXZ 
chain in Eq. (|2.2p with fii = /i (const.) is sketched in 
Fig. m At zero chemical potential, the chain is in its 
gapless, power-law correlated XY phase for —1 < 7 < 1. 
For 7 > 1 (7 < — 1), the spin chain assumes long- 
range Ising antiferromagnetic (ferromagnetic) order and 
the spectrum gaps out. We note that the thickness of 
the chemical potential window over which power-law XY 




FIG. 4. T = ground state phase diagram for the XXZ 
chain.— Red dashed lines correspond to constant fermion den- 
sity contours. 



order occurs (at fermion densities between and 1 per 
site) narrows to zero upon approaching the ferromagnetic 
transition at 7 = — 1. 

For the quantum quench studied here, the initial and 
final lattice Hamiltonians are given by 



H^'^ = - J 



(2.3a) 



E4 



Ci+l 



H.c. 



2MaE(-l)''5"i 



(2.3b) 



We assume periodic boundary conditions in a chain of 
N = L/a sites, with a the lattice spacing, so that cat+i = 
ci. We always take N to be an even integer. The initial 
Hamiltonian TJ''^ is tuned to reside in its gapless XY 
phase, so that —1 < 7 < 1. 

In Eq. ([Oa|) above. 



(0) 



QVtt ^(7) 
A ^(7)' 



(2.4) 



is the localized chemical potential used to introduce a 
particle density inhomogeneity near the center of the 
chain; we have introduced the spatial coordinate Xi = 
(i — N/2)a such that Xjv/2 = smd = o strad- 

dle the chain center. The parameters Q and A set 
the "strength" and width of the Gaussian potential, re- 
spectively. Two additional parameters that enter into 
Eq. (|2.4p are the sound velocity u and the Luttinger pa- 
rameter K. These coefficients completely determine the 
character of the low-energy field theory description of the 
XXZ chain in its critical XY phase, in equilibrium. In the 
absence of an external chemical potential, u and K can 
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be obtained from the Bethe ansatz, yielding^ 



IS 



u{'-f) — Ja 

Kin) 



arccos(7) ' 



2[7r — arccos(7)] ' 



(2.5a) 
(2.5b) 



such that u(0) = 2Ja = vp and K{Q) = 1. Here, vp de- 
notes the band Fermi velocity at half-fiUing with M = 0. 
We have included the ratio u/K \x\ the definition of the 
local potential so as to keep the initial induced density 
inhomogcncity approximately constant with varying in- 
teraction strength 7. 

The spectrum of H^f^ [Eq. ((23b)) ] is 



±2J 



\/cos2(fca) + (Maf. (2.6) 



The staggered potential, which doubles the unit cell, in- 
troduces a bandgap in the spectrum ed, kp = 7r/2a with 
magnitude Eg = AJMa. 

Based on the analysis of the continuum sine-Gordon 
quench in Ref. |43| . we expect the post-quench system re- 
sponse to be governed by the dynamical exponent. 



-(7) = \ 



1 



2[arcsin(7)]^ 
Ti'^ + 27r arcsin(7) 



- 1 



(2.7) 



Such an interaction-dependent exponent characterizes 
the (critical) power-law behavior exhibited by correlation 
functions in gapless ID quantum systems that possess a 
low-energy LL description. At 7 = (non-interacting 
quench), a assumes its minimum value of zero. At 7 = 1, 
on the precipice of the instability to Ising antiferromag- 
netism, a = 1/4. By contrast, a diverges upon approach- 
ing 7 = — 1 from above. 

In what follows, we set J = 1 and a = 1, thereby 
measuring energies in units of the transfer integral J and 
distances in units of the lattice spacing a. Our observable 
of interest will be the time-evolved density at each site 
of the lattice, i.e. 



p{t,Xi) = {5ni{t)) 



|e^«'"*4c,e-^'"*|0)-l/2. 



(2-8) 

where |0) is the ground state of the initial Hamiltonian 



+ [(-ir + (-l)^]G(2)(t,z-j) 
+ (-l)^+JG(3)(f,z-j)} 



JV 



= ^5,j(t)c,(0), 



(2.9) 



where Cj(0) denotes the Schrodinger picture operator, 
and 



N/2 



with 



TV 



N 



G^^\t, k) = cos{Ekt) - i-^ sm{Ekt), 
Ek 

o />/ 

G(2)(i,fc) = _j sin(i?fet), 

&^\t, k) - cos{Ekt) + sm{Ekt). 

Ek 



(2.10) 

(2.11a) 

(2.11b) 
(2.11c) 



In Eq. (|2.1ip . = — 2cos(fc) and Ek was defined by 
Eq. p.6p . The post-quench dynamics of the number den- 
sity are subsequently given by 

N 

p{t,x,)= [g*J{t)g^J'{t)c{x„x,,)], (2.12) 

where all information about the initial state is encoded 
in the static correlation function (single-particle density 
matrix) 



C(a;„x,v)s(0|c](0)c,.(0)|0); 



(2.13) 



|0) denotes the ground state of H^^K 

For the special case of the non-interacting quench, 
C{xj,Xji) is obtained by diagonalizing an N x N ma- 
trix. We denote the single particle Hamiltonian implied 
in Eq. (|2.3al) with 7 = by /i. In this case, the correlator 
is given by 



C{x.i,Xj')= UP{-hD)U^ 



3,0' 



(2.14) 



B. Dynamics; non-interacting quench 

For both the interacting (7 7^ 0) and non-interacting 
(7 = 0) quenches, the dynamics arc obtained by solving 
the Hcisenberg equation of motion for the annihilation 
operator ci (t) at site i using H'-f^ [Eq. ((2jb|) ]. Thercsuh 



where WhU = ho diagonalizes the single particle Hamil- 
tonian, and P{—ho) projects onto the (filled) nega- 
tive energy states of the diagonalized hp)- Combining 
Eqs. (|2.14p and (j2.12p gives the formal solution to the 
non-interacting quench. In practice, because of the in- 
homogeneity, we compute the single-particle matrix in 
Eq. (|2.14l) numerically. 
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III. CONTINUUM VS. LATTICE 



where the components '^p,=± (<, x) are defined via 



A. Warm-up: Relativistic wavepacket dynamics of 
a single massive Dirac particle 

Before turning to the continuum sine-Gordon quench, 
we pause to consider a toy problem: the time evolu- 
tion of a Gaussian wavepacket for a single, massive Dirac 
fcrmion in ID. This material is standard, but we include 
it to emphasize several important points regarding dis- 
parate regimes of relativistic wave propagation, and to 
clarify the similarities and differences between single par- 
ticle wave packet mechanics and the many particle quan- 
tum quench problem studied in this paper. 

The post-quench, band insulator Hamiltonian in 
Eq. (|2.3b[) exhibits a gap Eg = 2vpM, centered at 
k ~ kp = tt/2. Linearizing and truncating the band 
structure to modes near kp, one obtains 



where the single particle Hamiltonian is given by 



h = vp 



dx 



(3.1) 



(3.2) 



and the 2-component Dirac spinor ^/'(x) has the Fourier 
transform 



^(fc) 







e-*'^/4c(A: + kp) 






e*'^/4c(fc - kp) 



(3.3) 



with < < K <^ kp (A is a momentum cutoff). The 
components V'l and '02 denote right- and left-movers in 
the massless limit. In Eq. p.2|) . we have introduced a set 
of Pauli matrices {ct^''^''^} acting in the pseudospin space 
of tj}^ The Fermi velocity vp = 2; below we absorb it 
into the primed time, 



t' = Vpt. 



(3.4) 



In this section we take the system size L — > oo. 

We assume a Gaussian initial wavefunction for a par- 
ticle in its rest frame. 



(7rA2) 



2^l/4 



*o,i 

*0,2 



(3.5) 



with l^'o.iP = 1*0, 2p = 1/2. In this equation and the 
ones that follow, 4* denotes a single particle wavefunc- 
tion; its time evolution is determined by h in Eq. p.2p 
via the Schrodinger equation. It is useful to write the 
solution at times i > in two different ways. One way is 



*(t,x-) EE «'+(*, x) + *-(t,a;). 



(3.6) 



AA'/ 
(47r)3/4 



—i^A{t',x)cosh\z—iiZ(-\{t'^x)\ 



& 



(3.7) 



An alternative representation for ^'(t, x) is given by 



■^{t,x) 



1 



(7rA2)i/* 
+ y dye~ 



*o.2e"(-+*')'/2'^' 



-G<^'(t',y) G(3'(t',a) 



*0.2 

(3.8) 



where 



M 

T 

M 

T 

M 

T 



t + y 

Jo [A{t, y}] , 
t-y 



Ji[A{t,y)] 



(3.9a) 
(3.9b) 
(3.9c) 



denote the Af-dependent components of the Green's 
functions [from the continuum limit of Eq. (|2.1ip ]. In 
Eq. dXT]), tanh(zo) = x/t\ while A{t',y) = M^f^ - y^. 
Eqs. p.6p and p.7p follow from the momentum eigen- 
state expansion for the time evolution operator, while 
Eq. (j3.8l) obtains from the real space propagation am- 
phtude. In Eq. (|3.9p . the symbols J{o,i} denote Bessel 
functions of the first kind. 

A basic consequence of relativistic quantum field the- 
ory is that a single particle cannot be confined to a region 
smaller than its Compton wavelength 1/M. Localization 
to smaller scales induces particle energies in excess of the 
mass gap; in a many-particle theory, this typically leads 
to pair production out of the vacuum. 

In single particle relativistic wave mechanics, one in- 
stead finds qualitatively different behavior for initial con- 
finements A > 1/M ( "non-relativistic" ) and A 1/Af 
("relativistic"). We consider first the non-relativistic 
case. For MA ^ 1, the argument of the exponential 
in Eq. (|3.7p can be expanded to quadratic order in z. In 
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this approximation, one obtains 



where 



1 



40FA5(t') 



Sit') 
2f 



-cos [2Mt' - <j>Q{t',x)] 
-sm[2Mt' - (f>o{t',x)] 



_^l/[MA5{t')Y+2t'x/M^A*5^{t') 
1 



X < 



X 



1 + ^cos [2Aft'-(/)+(t',x)] 



_^l/[MA5{t')]^~2t'x/M^A'^S^{t') 
1 



X 



1 + ^cos [2Mt' ~4,.{t',x)] 
t' 



-sin [2Mt' - (j>_{t',x)] 



(3.10) 



where the scale factor S{t) = a/T+IVMW; 
(/>{0,+,-}(i7 a;) denote some phase factorsi^ For the ini- 
tial spinor components in Eq. p.Sp . we have made the 
choice ^o,{i,2} = exp(=F«7r/4)/-\/2, so that '^oi^) is h> 
variant under time-reversal and parity operationsi^ 

To the lowest order in 1/MA, Eq. (|3.10p reduces to 
the usual non-relativistic formula 



\^it,x)\' 



At smaller values of l/7\/A, the oscillatory character 
of Eq. p.lOp becomes important, and one observes the 
"Zitterbewegung" phenomenon: the evolving probabil- 
ity density exhibits an undulatory envelope that beats 
at the "interband" frequency w = 1/2M. These oscilla- 
tions occur because the eigenstate synthesis of the initial 
Gaussian [Eq. p.Sp ] requires larger contributions from 
negative energy states as the width A is narrowed. 

In the ultrarelativistic limit M = 0, Eq. p.Sp implies 
that 



\^it,x)\' 



1 



2V^A 



,^{x^t'f/A^_^^-{x+t'f/ 



A^ 



By contrast, when < MA ^ 1, the propagation is rela- 
tivistic but dispersive. In the long time limit, the second 
term on the right-hand side of Eq. (|3.8p is dominated by 
the diagonal Green's functions. Close to the right light- 
cone edge |a; — t| < A, ior t' ^ t'^-^^ one obtains 



1 



-{x-t'f/2A^ 



+ dispersive background 



/ t' x-t'\ 










(3.11) 



t' = l/2M^A. 



The function 



dw 



1 



Ji{w) 



vanishes in the limit a — !■ oo. 

In the XXZ chain quantum quench studied in this pa- 
per, it will prove essential to distinguish rclativistic vs. 
non-relativistic initial conditions using the width A of the 
Gaussian chemical potential inhomogeneity in Eq. (|2.4p 
and the bandgap parameter M in the post-quench Hamil- 
tonian Hf [Eq. (j2.3bp ]. Examples of single particle 
non-relativistic and rclativistic propagation are shown in 

Fig.ia 



B. Sine-Gordon quench and "supersolitons" 

We now consider the continuum limit of the XXZ 
quench defined by H^"^ and H^^^ [Eqs. (f23a|) andJpEl). 
This problem was previously analyzed in Ref. l43l In 
this section, we provide the solution to the sine-Gordon 
quench and a brief recapitulation of the results found in 
Ref. 0. In Sec. IIII D[ we consider the modification of 
these results due to the presence of irrelevant operators 



a) 




FIG. 5. Examples of single particle, massive Dirac equation 
wave packet propagation, obtained via numerical integration 
of Eqs. (|3.6p and p.7p . The case (a) corresponds to a "non- 
relativistic" initial condition, the Gaussian in Eq. (|3.5fl with 
A = 5/M. The "rclativistic" case is illustrated in (b), with 
an initial A = 0.5/Af. Here M = 1, and data is shown at 
times t = 0, 15, 35, and 55; fainter (bolder) traces depict 
earlier (later) times. 



9 



(i.e., lattice-scale details left out of the renormalizable 
continuum field theory). 

The massive Dirac continuum limit for the final 
state Hamiltonian H^^'^ was derived in the last sec- 
tion, Eqs. (|3.ip - (|3.3p . Since this Hamiltonian is non- 
interacting, we can construct a formal solution to the 
quench dynamics by solving the Hcisenberg equations of 
motion for the Dirac spinor tp(t,x) [c.f. Eq. (|2.9p ]. The 
result is 



dyJ{t\y)i,{Q,x-y), 



where 



J{t,y) = Mt,y)+JM{t,y) 
5{y - t) 



Mt,y) 

Juit^y) 



5{y + t) 
&'Ht,y) &^Ht,y) 



(3.12a) 

(3.12b) 
(3.12c) 

(3.12d) 



Eq. p.l2p is identical to the propagation amplitude tran- 
scribed in the last section [Eq. p.Sp ]. after replacing the 
single particle wave function vE'o(x) with the Schrodinger 
picture field operator ■0(0, cc). In these equations, t' — 
vpt [Eq. (inH)], while the Green's functions (^(^■^'^'(t, y) 
were defined by Eq. p.9p . 

The post-quench fermion density is given by 



p{t, x) 



dyi 



where all information about the initial state is encoded 
in the correlation function 

C\{xi,X2) = (O[V'"^(O,xi)0,(O,X2)|O). (3.14) 

In Eqs. (|3TT3l) and (f3J4|) . the indices r,s e {1,2}; re- 
peated indices arc summed. 

The pre-quench system is described by the ket [0), 
which is taken as the ground state of the Luttinger liquid 
Hamiltonian 



dx 



dx 



,(0) 



(x) : iI^'^iIj : +2^vf '■ i^'^ tp '■ 



(3.15) 

Eq. (|3.15p gives the continuum limit of i/^'^ in Eq. p.3ap . 
after discarding all irrelevant operators; here, fi^^^x) rep- 
resents the long wavelength, continuum approximation 
to the lattice potential /xl"-* The symbol : . . . : denotes 
normal-ordering. 

Using abelian bosonization rulesj^i^ we rewrite 
Eq. ((XT5|) as 



dx 



uK 



dx 



+ 



2K 



dx J 



^(")(x) d0 
y/n dx 
(3.16) 



In our conventions, the fermion current components are 
bosonized as 

(3.17) 

and satisfy [J°{x), J^{x')] = -{i/Tr){d/dx)S{x - x'). The 
sound velocity u and the Luttinger parameter K in 
Eq. (|3.16p are given by 



1 



1 



47 



(3.18) 



With vp ~ 2, these agree with the Bethe ansatz results 
in Eq. (|2.5p only to the first order in 7. 

When expressed in terms of the boson variables, the 
post-quench massive Dirac Hamiltonian in Eq. (|3.ip be- 
comes the sine-Gordon model 



H^f^ = VF I dx 



1 

2 \dx 
M 



1 I'M 

2 [d^ 



— cos ( VAttO 



dy2 J\t',y^)j{t',y2) 



X C^ix - yi,x - y2), (3.13) [: 



(3.19) 



The variable a appearing in the prefactor of the cosine 
term carries units of length, and is formally introduced 
by the bosonization procedure.—"^ 

While -ff'*' [Eq. (|3.16p ] assumes a non-interacting 
form when expressed in boson variables, H^^^ becomes 
the non-linear sine-Gordon theory. By contrast, H^^^ 
Eq. p.ip ] is non-interacting in terms of the Fermi field 
-0, while /f^*) [Eq. p.lSp ] incorporates four fermion in- 
teractions. For a quench with 7 7^ 0, there is no common 
language in which both ij(') and H'^f^ can be simulta- 
neously expressed as non-interacting Hamiltonians. We 
refer to this generic scenario as the "interacting" quench 
in the sine- Gordon theory. We reserve the appellation 
"non-interacting" for the exceptional case with 7 = 0, 
where both and TJ^-^^ are bilinear in fermions. The 
post-quench dynamics exhibited for each case are differ- 
ent, as discussed below. 

Expressing the fermions in Eq. p.l4p as vertex opera- 
tors in the bosonic languagcji the initial state correlation 
function components evaluate to 



C\{xi,X2) = 



ICNd 



■ exp 



27r 

sgn(a;i - X2) 



,£i / dyf,^°\y) 



\Xl - X2 



<T + 1 



(3.20a) 



C 2(2^1, 2:2) = — ^ — exp 



sgn(a;i - X2) 



\Xi - X2 



<T + 1 



C\{xi,X2) =C\{xi,X2) = (L ^ 00). 



(3.20b) 
(3.20c) 



In these equations, the external chemical potential man- 
ifests in a gauge "string" due to the axial anomalyj^i^ 
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The coefficient a was introduced in Eq. p.l9|) : the pa- 
rameter c^v is a numerical normalization constant*^ The 
off-diagonal components of g{xi,X2) {r ^ s) vanish in 
the thermodynamic limit L — > oo. 

The essential character of the initial Luttinger liquid 
state is encoded in the dynamic exponent cr, defined as 



(3.21) 



For the non-interacting quench, K ^ 1 and 7 = tr = 0. 
By contrast, any ^ 1 (7 7^ 0) gives a > 0,^ Eq. ([5:^ 
implies that a is twice the anomalous scaling dimension 
of in the initial LL ground state. 

Using Eq. p.20p . the integrals appearing in the fi- 
nal expression for the post-quench density expectation 
in Eq. p.l3p are ultraviolet convergent for < cr < 1. 
Over this range of initial conditions, we obtain a cutoff- 
independent prediction for the post-quench evolution of 
the number density in the continuum sine Gordon field 
theory. Similar expressions with identical convergence 
properties may be obtained for the kinetic and potential 
energy densities due to the inhomogeneous initial state 
chemical potential4^ 

The general characteristics of the long-time density dy- 
namics implied by Eq. (|3.13p and (|3.20p for a generic 
initial state ^^^-^ {x) were discussed in Ref. In this pa- 
per, we restrict our attention to the waves induced by a 
localized, Gaussian initial inhomogeneity. 



-ni^Hx) 



(3.22) 



Combining Eqs. (|3^ . (|3^ . and ((3T3)) gives an exact 
integral expression for the post-quench density expecta- 
tion evolution after the sine-Gordon quench. In the long 
time limit t 3> l/vpM, the requisite integrals yield to 
a systematic asymptotic analysis, as explained in Ap- 
pendix El One thereby obtains the exact leading asymp- 
totic behavior 



p(t,x) 



Q 



'20FA 

_ Q r(i-g) 
2Ar(i±^) 

+ {x ^ — x}. 



V2A 



a/2 



X-t' 



(3.23) 



where 



F,{z) = cxp(-zV2)Z?,/2(V2z), (3.24) 



and Di,{x) denotes the parabolic cylinder function. In 
Eq. p.23|) . we have used the explicit expression for the 
normalization constant Regardless, for cr > (in- 

teracting quench) the prefactor of the second term in 
Eq. p.23p is in some sense arbitrary, due to the a factor. 
This ambiguity can be resolved if a conventional normal- 
ization is adopted for the vertex function correlators in 




FIG. 6. The right-moving "supersoliton" obtained in the 
interacting sine-Gordon quench. The number density evolu- 
tion after a Luttinger liquid to insulator quench is depicted 
for a Gaussian initial density profile (heavy black line), with 
(T = 0.7, A = 3, and M = 15/16, obtained via numerical 
integration of the exact bosonization result [Eqs. (|3.13l) and 
(|3.20p . using Eq. p.22[l ]. Time series for two different Q are 
plotted; the densities are normalized relative to these. The 
evolution is reflection symmetric about x = Q. 



Eq. The derivation of Eq. ((X^ is sketched in 

Appendix \K\ 

For the interacting quench with < cr < 1, Eq. p.23p 
describes the propagation of right and left moving "su- 
persolitons" launched from the Gaussian initial condi- 
tion, in the long time limit. A right-moving supersoliton 
is depicted in Fig. [SI From the equation, it is evident that 
the supersoliton does not disperse. In the long-time limit, 
the response to the initial chemical potential (and thus 
the initial state density inhomogeneity) is linear, regard- 
less of the strength of Q in Eq. (|3.22p . The supersoliton 
features an amplitude that grows in time as power law, 
with growth exponent cr/2; subleading terms neglected 
in Eq. p.23|) decay for cr > 0. The peculiar shape of 
the supersoliton implied by Eq. p.24p obtains because 
the quench kernel effectively takes a fractional deriva- 
tive [d/dxY^"^ of the input profile.— The total number 
fluctuation induced by the inhomogeneity is conserved 
by the supersoliton, since the second term in Eq. p.23p 



integrates to zero over x G 



59 



Finally, we note that 



Eq. (|3.23p holds for generic MA: the supersoliton arises 
for both "relativistic" (MA ^1) and "non-relativistic" 
(MA < 1) initial density profiles (c.f. Sec UlT^ . 

By contrast, the non-interacting quench with a = 
exhibits no amplification. For MA ^ 1 ("non- 
relativistic" initial condition), one finds simple disper- 
sive broadening, qualitatively similar to the single parti- 
cle wavepacket spreading in Fig. [SKa). Examples of non- 
interacting quenches with non-relativistic initial condi- 
tions are shown in Fig.[7l For the non- interacting quench, 
the response is given entirely by terms neglected as sub- 
leading (for cr > 0) in Eq. (|3.23p : indeed, the right-hand 
side of this expression vanishes for cr = 0. 

Non-interacting quenches with "relativistic" (Af A ^ 
1) initial conditions exhibit a different behavior, qualita- 
tively similar to the single particle wavepacket evolution 
depicted in Fig. E^b): the initial Gaussian density bump 
blows apart into left and right-moving wave trains, with 
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leading edges that rip along lightcone. In this sense, the 
non-interacting quench with Af A <^ 1 behaves similar to 
the supersoliton, which also propagates relativistically. 
At short timescales t < l/vpM, the interacting and non- 
interacting quenches in fact exhibit qualitatively simi- 
lar dynamics for relativistic initial conditions. However, 
the non-interacting quench evolution shows no amplifica- 
tion in the long time limit, and the generated wave train 
exhibits no static, non-dispersing structure. In the nu- 
merical results for the XXZ chain quench presented in 
Sec. IIVI we will restrict our attention to non-relativistic 
initial conditions, in order to avoid possibly confusing the 
supersoliton with the trivial (and essentially single parti- 
cle) effect of squeezing the initial density wavepacket to 
a width narrower than the Compton wavelength. 

Finally, we note that setting M = in Eq. p.23|) gives 
the result appropriate to a Luttinger liquid to Luttinger 
liquid quench-the initial density disturbance is merely 
propagated along the light cone without dispersion, as 
expected for dynamics generated by a critical state. This 
case was previously considered in Ref. H^. 

The physics of the interacting and non-interacting 
quenches are fundamentally distinguished by the advent 
of quasiparticle fractionalization in the interacting case. 
The interacting nature of the pre-quench initial ground 
state |0) relative to the post-quench Hamiltonian ij'-^-* is 
implied by Eq. p.lSp with 7 7^ 0, which expresses -ff''' 
in terms of the "final state" fermion V'- The presence 
of interparticle interactions in iJ^'-' means the fermion ip 
is not a "natural" propagating degree of freedom in the 
initial state Luttinger liquidi^ Equivalently, the "quasi- 
particles" of the LL carry a fraction \/K of the ip fermion 
number charge; we say that the initial LL state is frac- 
tionalized with respect to the final band insulating state. 
This notion is made explicit in Appendix |B] Fractional- 
ization due to the presence of interparticle interactions 
is a ubiquitous feature in ID, responsible e.g. for spin- 
charge separation in quantum wires.— 




FIG. 7. The number density evolution as in Fig. [6l but for 
the non-interacting quench K = 1 {a = 0). The initial bump 
(heavy black line) has area Q = 0.1 in the main figure and 
Q = 1 in the inset; in both cases A = 3 and Af — 15/16. The 
evolution is reflection symmetric about a; = 0. Now there is no 
fractionalization of the initial LL quasiparticles with respect 
to the insulator and, consequently, the dynamics are simply 
dispersive with no supersolitons or inhomogeneity growth. 



We interpret the supersoliton and the amplification ef- 
fect that arises in the sine-Gordon quench for the inter- 
acting case as a spectroscopy of the initial LL state.— 
The key ingredients are quasiparticle fractionalization of 
the initial state relative to the excitation spectrum of 
the post-quench Hamiltonian, and fact that the post- 
quench Hamiltonian is gapped. In the next section and 
in Appendix \C\ we show that fractionalization leads to 
an anomalous momentum dependence in the Wigner dis- 
tribution function of the excited post-quench quasiparti- 
cles, due to the inhomogeneity. The low-energy disper- 
sion of the gapped final state characterized by Af trans- 
lates this into a divergent velocity distribution, giving rise 
to the supersoliton. By contrast, for a non-relativistic 
density profile with Af A ^ 1, Pauli-blocking suppresses 
the excitation of large velocities in the non-interacting 
quench. The distinction arises due to the long-distance 
behavior of correlations in the initial state, and is not 
destroyed by lattice effects. At the same time, we will 
see in Sec. IIII D] that the advent of the lattice does mod- 
ify the post-quench dynamics, but in way that can be 
parametrically controlled by the system size. 



Quasiparticle distribution functions: continuum 
and lattice quenches 



We consider the post-quench distribution of final state 
quasiparticles, in the sine-Gordon and lattice quenches. 
Time evolution is generated by the final state Hamil- 
tonian, which is translationally invariant and non- 
interacting in terms of band fermions for both the lat- 
tice [ijW), Eq. dlSbl] and continuum [H'^^\ Eq. ((3Tl) ] 
theories. The global momentum distribution of excited 
quasiparticles induced by the quench in each case consti- 
tutes a static quantity, which does not encode informa- 
tion about the density inhomogeneity. The physics of the 
supersoliton resides in the local Wigner function, which 
is discussed subsequently. 



1. Global distribution function 

We consider first the continuum sine-Gordon quench. 
The final state Hamiltonian H^^'> in Eq. (|3.ip and 
can be rewritten as 



H 



(/) 



dk 
2^ 



(3.25) 



where (6^) annihilates a particle (hole) with momen- 
tum fc, and 



£k 



(3.26) 
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The particle and hole operators arc related to ip via 



ak 



bL. 



1 

is{k) 



— is(— fc) 



(3.27) 



where 



sik) 



vfM 

We define the occupation numbers 

n+{k) = (0|4afc|0) 



n-{k) 



- /ni;,t 



blhk\0), 



(3.28) 



in which |0) denotes the ground state of the pre-quench 
Hamihonian H'^^\ Eqs. ((XT5|) and (PH)) . 

For the translationally invariant case //'■'^•' [x) = 0, one 
can show that 



n+{k) 



where 



F{k) = (0|V'^(fc)Vi(fc)|0) 



(3.29) 



(3.30) 



and F{-k) = -F{k). The form of Eqs. (jX^ and (P^O)) 
follows from Eq. (|3.27p and the imposition of the sum rule 
(canonical anticommutation relations) upon correlators 
of the fermion components ipi{k). 

In the case of the non- interacting quench (cr = 0), one 
finds J'(fc) = -(l/2)sgn(A:), so that 



n±{k) 



1 



+ A/2 



(3.31) 



This occupancy factor peaks to a value of one-half at 
fc = 0, and decays as M^/fc^ for |fc| > M. The density of 
particles or holes excited by the quench is thus ultraviolet 
finite and equal to M/2n. The associated kinetic energy 
density is given by the difference of the pre- and post- 
quench Hamiltonian zero point energy densities, and is 
logarithmically divergent. 

The calculation for the interacting case is more subtle, 
due to an ultraviolet divergence. We must compute 



(0|^it(fc)Vi(fc)|0) 



cnol 
27r 



dx 



e [a -t- ix) 
Ja^^x^f^' 



(3.32) 

which is the Fourier transform of the initial state cor- 
relator in Eq. p.20a[) with = 0, retaining the soft 
cutoff aX The prefactor cat appears explicitly in Ref. [sgI . 
Performing an expansion in ka and extracting T{k), we 
finally obtain 



\k\ r(V) 

Vfc2 + Af2r(i±2:) 



\k\a 



0(\k\a), 



(3.33) 




+ 



0.5 
0.4 
0.3 
0.2 
0.1 
0.0, 



JL ' 1 l' 1 ' 


1 ' 1 ' 

(c) : 
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FIG. 8. Occupancy of momentum states in the upper (con- 
duction) band after the lattice XXZ quench (obtained by 
DMRG), for the four values of the dynamical exponent a 
given in the legend. The occupancies are plotted for four 
different choices of the band-gap parameter, M = 3/40 (a), 
1/8 (b), 1/4 (c), and 3/8 (d). In each subplot, the dashed 
vertical line marks the wavenumber fcmax(Af) at which the 
final state Hamiltonian band group velocity is maximized, 
u(fcmax;M) = max[d£fe/dfc] = Wmax(M); see Sec. IIVB2I 
Eqs. (113} and gH]). 



valid for < cr < 1. Eq. p.33|) holds only for \k\a small, 
where the second term on the right-hand side trails the 
first. At such wavevectors, n± is enhanced relative to the 
non- interacting case in Eq. p.3ip . indicating that the 
interacting quench induces a stronger excitation of the 
post-quench quasiparticles. Clearly Eq. p.33p becomes 
unphysical for sufficiently large |fc|; the global distribu- 
tion function cannot be uniquely defined (i.e., its value 
will depend upon the regularization procedure) in the 
continuum, interacting sine-Gordon quench. 

In Fig. [8l we exhibit n+{k) for the lattice quench in 
a finite size system of 202 sites, obtained via numer- 
ical density matrix renormalization group calculations 
(see Sec. IIVI for details). The occupancy is defined as 
in Eq. (|3.28p . except that the continuum state |0) is re- 
placed by |0), the ground state of H^^^ in Eq. (|2.3a[) : ak 
now denotes the lattice conduction band annihilation op- 
erator. 



Wigner function, fractionalization, and the origin of the 
supersoliton 



Eqs. p.3ip and p.33p have been calculated for the 
case of the homogeneous quench. In the infinite sys- 
tem size limit, these equations also apply in the presence 
of an arbitrary initial state chemical potential /x^''^(x) 
that vanishes faster than 1/x as |a;| — >■ cjo. The effects 
of the inhomogeneous density profile in the initial state 
can be extracted from the "local" (Wigner) distribution 
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functioni^ 

The main idea is conveyed by the ground state Wigner 
function for the right-moving fermion ipi{x) = TZ{x) in 
the inhomogeneous Luttinger hquid, defined as 



nnik] R) 



; J dxd e-^'^""" (0|7^^ {R-^)n{R+ ^) |0) 
ge»'«(0|7tt (fc_|)7l(fc+f)|0). 



(3.34) 

Here |0) is the ground state of 5^*' in Eq. (I3.15p . which 
has the density expectation lfi\R^ {x)TZ{x)\Q) — pq{x)/2, 
where 



fi^'»{x). 



(3.35) 



The global momentum distribution and the position 
space density expectation value can both be extracted 
from the Wigner function (Appendix [C|) . Although we 
employ it here to gain intuition about the local momen- 
tum profile induced by Pq{x), strictly speaking n-]^{k] R) 
cannot be interpreted as a probability distribution^ (it 
can take negative values), because momentum and posi- 
tion are canonically conjugate quantum observables. 

We let Sn-ji{k;R) denote the linear response due 
to po(x), after subtracting off the homogeneous 
(global) background. Using the correlation function in 
Eq. p.20a|) . we find 



5nu{k] R) 



2tt q 



sgn [k + - j 

- sgn {k - 1^ Ga 



k + 



k-1 



■A 



(3.36) 



where 



dy 



y (y2 + (^2) 



The parameter C is an ultraviolet regularization length, 
introduced here for later use in the context of the lattice 
quench defined in Sec. Ill Al (see Sec. IIIIPp . The pure 
continuum theory has C = 0, for which Eq. p.37p is con- 
vergent when < cr < 1. 

For the special case of a non-interacting Fermi gas {a = 
0), ^odphC) ^ ^/2, independent of so that 
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ign (fc + I) - sgn (k - I) 



(3.38) 

This expression vanishes for \q\ < 2\k\, i.e. unless the 
creation and annihilation operators in Eq. (|3.34l) carry 



momentum with opposite signs. This is a simple conse- 
quence of Pauli blocking. The result can be understood 
via perturbation theory: let |0) = |0)o 4- \S); |0)o denotes 
the homogeneous vacuum, while \5) gives the response 
to po{x). To first order in \S), Eq. p.34p vanishes unless 
the product W {k ~ f) ^ f) creates a particle-hole 
pair in |())o (acting to the left or to the right); Eq. p.38p 
obtains from the overlap between this state and | S) . 

We consider a po{x) localized in position space, of char- 
acteristic width A. For \k\ > 1/A, Eq. implies 
that the fc-dependence of the Wigner function is slaved 
to follow that of po{q), with q ^ 2k. For a Gaussian inho- 
mogcneity, this means a Gaussian fall off of the Wigner 
function in k. No matter how wide or narrow the initial 
packet is made, the large- fc asymptotic is always strongly 
suppressed. 

By contrast, the situation for the interacting Luttinger 
liquid is quite different, Eq. p.36|) with cr > 0. Then, the 
kernel Ga-i\p\',C) depends upon \p\, and allows a contri- 
bution to 6n-jz{k;R) from \q\ < 2\k\, violating the Pauh 
blocking condition. For \k\ ^ 1/A, the dominant contri- 
bution comes from \q\ <^ 2\k\, and the (/-integration gives 
Sn-]z{k; R) oc p{R). The fc-dependence comes entirely 
from dC/o-dfcl; C)/o'|fc|i ^^^d is independent of the initial 
inhomogencity profile: the R- and fc-dependencies of the 
Wigner function factorize. For < a < I and C = 0, the 
leading term in the ^ 1/A limit of Ec^. p.36p goes as 

5nn{k;R) - c(cr) a'"|fc|'""Vo(-R), (3.39) 

where the prefactor satisfies c(0) = 0. 

In the sine-Gordon quench, the initial "momentum 
distribution" implied by Eq. p.36p or (|3.39p translates 
into a "velocity distribution" through the massive post- 
quench dispersion relation in Eq. p.26p : details are pre- 
sented in Appendix [Cl In the non-interacting case, for 
a non-relativistic initial condition with MA ^ 1, only 
small velocities v < vp/MA are excited [Eq. (|C7|) in 
Appendix [C] . By contrast, an interacting quench with 
Eq. (|3.39p induces a non-integrable divergence in the 
velocity distribution at the "speed of light" v ^ vp 
[Eq. (jCSap j. signaling the presence of the non-dispersing 
supcrsoliton (recall that velocity is conserved by the post- 
quench Hamiltonian). Thus, the supersoliton arises due 
to the particular way in which quasiparticle fractional- 
ization evades Pauli blocking in the initial LL ground 
state. 

To gain further insight, consider a many-particle prod- 
uct state in a relativistic quantum theory, e.g. N parti- 
cles are placed into N plane wave states with momenta 
{ki}. In the thermodynamic limit N — !■ oo. the system 
is described by a continuous distribution function n{k), 
a well-defined (classical) observable for a product state. 
Suppose further that the corresponding velocity distri- 
bution n(v) exhibits a delta-function-like singularity at 
V = vp- Then, a fraction of the density (determined by 
the weight of the singularity) at each point in space is 
translated at the speed of light. In particular, a frac- 
tion of any initial density inhomogencity will propagate 
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at V = vf without dispersion. 

Even though this intuitive velocity distribution picture 
helps reveal the physical origin of the supersoliton, it is 
not entirely satisfactory. For example, a finite number 
of massive particles traveling at the speed of light im- 
plies an infinite kinetic energy. By contrast, although 
it propagates ultrarelativistically, the supersoliton car- 
ries finite total energy beyond that induced by the ho- 
mogeneous quench4^ The initial (inhomogeneous) Lut- 
tinger liquid is very far from a product state of the post- 
quench spectrum; indeed, the appearance of the expo- 
nent a in Eqs. p.33p and p.39p indicates that quantum 
coherence (entanglement) plays a dominant role in the 
fractionalized density dynamics of the quench. Further, 
the Wigner function Sn{k;R) is not really a phase space 
distribution function^ and the non-integrable velocity 
singularity in Eq. (|C8aP , Appendix [C] does not imply an 
extensive mass or energy flow. Instead, we interpret this 
divergence as signaling the supersoliton, an emergent, 
collective excitation of the post-quench non-equilibrium 
state that travels with velocity vp. 

Let us also note that the Wigner distribution post- 
quench is not a static object; the spatially-varying mo- 
mentum profile in Eq. p.39|) implies that the shape of 
the density wave can evolve. In the unregularized sine- 
Gordon quench, this is the amplification effect exhibited 
in Fig. m and Eq. ([3221) • 

The ultraviolet effects induced by the presence of a lat- 
tice cannot alter the fundamental distinction between in- 
teracting and non-interacting quenches, because the long- 
distance behavior of {0\TZ^ {0)TZ{x)\0) determines the effi- 
cacy of Pauli blocking in Eq. p.36[) through Eq. (|3.37p . In 
the next section, we will nevertheless see that the modifi- 
cation of the short-distance structure of the theory [e.g., 
C > in Eq. p.37p ] does influence the post-quench dy- 
namics. 

Further details about the post-quench Wigner function 
can be found in Appendix [Cj where explicit formulae are 
given for the associated velocity distributions in the non- 
interacting and interacting quenches, incorporating the 
effects of ultraviolet regularization. 



D. Irrelevant operators and UV regularization 

So far, we have focused primarily on the Luttinger 
liquid to band insulator quench in the continuum sine- 
Gordon model. For a Gaussian initial density bump, 
the leading asymptotic result for the long-time limit 
[Eq. p.23p ] predicts the emergence of the supersoliton: a 
non-dispersive, relativistically-propagating density wave 
with an amplitude that grows as t"'/^. This result ap- 
plies to the integrable sine- Gordon field theory, in the 
absence of additional perturbations. We have considered 
only a particular case by assuming the non-interacting 
post-quench Hamiltonian in Eq. p.ip . This corresponds 
to the special Luther-Emery (LE) point in the sine- 
Gordon phase diagramjii^i^ Away from this point, H'^^'> 



would acquire a four-fcrmion interaction as in Eq. (|3.15p ; 
bosonization links the sine-Gordon and massive Thirring 
models in the general casci^ We postpone a discussion 
of the effects of interparticle collisions in the post-quench 
evolution until the end of Sec. |Vl 

The sine-Gordon field theory arises as the low-energy 
description of many ID solid state and cold atomic 
systemsj^i^ including the XXZ chain quench intro- 
duced in Sec. EX] between and iJ^-^' in Eq. (lO]) . 
Details of the original "microscopic" formulation are ex- 
pected to appear in the effective low-energy field theory 
as irrelevant opcratorsi^ Irrelevant operators typically 
exert a negligible effect upon low-energy, long- wavelength 
properties in a zero temperature field theory. Finiteness 
of correlation functions (up to logarithmic divergences 
subsumed by renormalization) and insensitivity to irrel- 
evant operators go hand-in-hand.— 

By contrast, the influence of irrelevant operators upon 
the strong non-equilibrium dynamics generated by a sud- 
den quench remains largely unexplored territory. The in- 
corporation of generic perturbations destroys some spe- 
cial properties that may be enjoyed by a given renormal- 
izable theory, such as conformal invariance or, in the case 
of the ID sine-Gordon model, integrability. On general 
grounds, a non-integrable many-body system prepared 
in an initial, non-thermal state is expected to thermalize 
(presumably due to quantum chaotic dynamics) in the 
long time limit. 

In this paper, we do not mount a broad assault upon 
the important topics of intcgrability-breaking perturba- 
tions and thermalization. Even in equilibrium, the im- 
pact of irrelevant operators and integrability on corre- 
lation functions at non-zero temperature remains a con- 
tentious issue4^i^"— Here, we limit our focus to the post- 
quench wave train dynamics exemplified by the super- 
soliton. In particular, we would like to understand how 
irrelevant operators, or cquivalently, lattice scale details 
and the presence of a finite ultraviolet cutoff, modify or 
suppress the supersoliton. Our considerations in this sec- 
tion will be used to interpret the numerical results for the 
XXZ chain quench presented in Sec. IIVI 



1. Irrelevant operators: some examples 

The XXZ chain quench introduced in Sec. Ill Al takes 
the ground state [0) of the XY phase Hamiltonian iJ^*^ 
in Eq. (|2.3a|) . and evolves this state forward in time using 
the gapped band insulator Hamiltonian TJ'-^^ defined via 
Eq. (j2.3b[) . In the continuum field theory limit, lattice 
microscopies induce the addition of irrelevant operators 
to the sine-Gordon model Hamiltonians if'*^ and H^^'^ 
[Eqs. ([XT5)) . ((XTB)) . and ([SIT]), ([XTO)) ]. We now enumerate 
a few examples. 

The least irrelevant operators {Oi{x)} invariant un- 
der continuum versions of all lattice symmetries (time- 
reversal, parity, lattice translational invariance) carry 
the scaling dimension Xi = 4 when added to the non- 
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interacting Dirac Hamiltonian in either Eq. p.l5|) (with 
7 = 0) or p.ip . We consider first the uniklapp interac- 
tion operator— 



Ou{x) 



1 



(7ra) 



■ cos 



2V47r6' 



(3.40) 



This operator appears as a lattice-induced modification 
(via iJ^*)) of the initial Luttinger liquid Hamiltonian 
-ff , Eq. ([XT5)) or ^ The dimension of the umklapp 

operator is xi = AK^ so that the associated coupling con- 
stant has dimension yi = 2(1 — 2K), where K denotes the 
Bethe ansatz Luttinger parameter in Eq. (j2.5bp i^ Thus 
the umklapp operator has yi{K = 1) = —2 at the free 
fermion point, while yi{K = 1/2) = at the threshold of 
the Ising antiferromagnetic order [Eq. (|2.5bp with 7 — > 1] . 
In our lattice quenches, we will focus upon 7 < 0, so that 
K > 1 and umklapps are strongly irrelevant. 

As a second example, we consider the effect of band 
curvature (at half- filling) , which gives the operator 



"3 \ , 



1 
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dx"^ 



di) 



d<t) 

dx 



de 

dx 



(3.41) 



©3 arises as a modification of both if^*^ and H'^^\ due 
to the cosine dispersion of the lattice model. 

We note that while the band curvature operator in 
Eq. p.4ip is bilinear in terms of fermions, both the umk- 
lapp and band curvature operators induce interparticle 
interactions in the boson language. This complication 
makes it difficult to determine the influence of either 
upon the interacting LL initial state |0). 

One can in principle treat non-bilinear irrelevant op- 
erators perturbatively, but several difficulties arise in at- 
tempting to account for their effects. First, the pertur- 
bation theory is badly ultraviolet divergent, and depends 
upon the way in which these divergences are regularized. 
A second, more serious (but intimately related) problem 
arises because the effects of irrelevant operators become 
strong at short distances. In the context of the quench, 
the goal is to construct the initial ground state correla- 
tor in Eq. (|3.14p , accounting for the effects of lattice scale 
details. These details should translate into a modifica- 
tion of the ideal Luttinger liquid correlation functions in 
Eq. p.20p at short distances. However, the effects of ir- 
relevant operators become strong in precisely this limit; 
the result is that perturbation theory breaks down, and 
a systematic accounting is only possible via an exact or 
approximate non-perturbative resummation. 



Nevertheless, we show that a regularized version of the 
continuum sine-Gordon quench can be constructed which 
gives a reasonably good match to our finite system size 
numerics presented in Sec. IIVI To motivate the regu- 
larization scheme that we employ, we consider the effect 
of a finite-ranged density-density interactiouii Instead of 
Eq. (|3.16l) . one has the Hamiltonian 



1 

+ 2^ 



dx 



dx dx' v(x 



d£Y 

dx J 



dx dx 



(3.42) 



where v{x) = (27/C) exp(— |a;|/^) is a Yukawa-type po- 
tential that integrates to 47, regardless of the range C,- 
The limit C ^ gives the purely local interaction imple- 
mented in Eq. (|3.16p . 

Unlike the non-bilinear boson operators associated 
with umklapp and band curvature effects discussed 
above, the finite-range interaction in Eq. p.42p can be 
treated non-pcrturbatively. The result is a modification 
of the Luttinger liquid correlator in Eq. (|3.20p : in the 
homogeneous limit with ^'■'^■'(a;) = 0, one obtains 



c\(,,o)=^^e:(:r). 



27r 



(3.43) 



where 



<L{x) 



- exp 



dq 



<j{q) [cos (qx) - 1] 



C+I3x^ 



7(0)/2 



In this equation. 



K{q) EE 



[Kiq)+K-\q)]-l, 
1/2 



1 



47 



1 



TT 1 + (Cq)2 



(3.44) 

(3.45a) 
(3.45b) 



The variable /3 in Eq. p.44p is some numerical constant. 
The effect of a finite interaction range C > is to re- 
duce the short-range scaling behavior {\x\ < C) of the ip 
fermion LL correlation functions in Eq. p.20p to that of 
free fermions. 



2. Regularized sine-Gordon theory 

A systematic approach to incorporating lattice scale 
details into the sine-Gordon quench would require the 
inclusion of all irrelevant operators with a given scaling 
dimension, say. This task is made difficult by the inter- 
acting nature of most such operators. The problem is 
compounded by the fact that the influence of all irrele- 
vant operators becomes strong in the ultraviolet, which 
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FIG. 9. The ("-regularized supersoliton obtained by numeri- 
cally integrating Eq. (|3.47p . Here we have set = 1, A = 6, 
and plotted data for four values of a. The interaction expo- 
nent a = 0.7. We have assigned M = 3/2A, so that the dy- 
namics reside within the "non-relativistic" transport regime, 
as discussed in Sec. HIT Al For all but the smallest value a 
depicted, the characteristic "s" shape of the supersoliton is 
identified at sufficiently long times. In comparison to the pure 
sine-Gordon model result shown in Fig. [6l the amplitude of 
the regularized supersoliton saturates at times t > ti^. Two 
inequivalent definitions for are provided by Eqs. (|3.49[) and 



is precisely the regime where lattice scale effects are ex- 
pected to manifest F^i 



In the following, we sidestep these difficulties with a 
phenomenological approach. Compare the lattice and 
continuum initial state correlation functions C{xj,Xj') 
and C'',(a;i,a;2) defined by Eqs. ([^1^ and i^J^. While 
the long-distance behaviors of these functions should be 
compatible, the short-distance behaviors clearly differ. 
The continuum LL correlation functions in Eq. p.20p 
exhibit a power-law divergence as xi — > X2 governed 
by twice the scaling dimension {a + l)/2; by contrast, 
the lattice correlator satisfies C{xj,Xj) = (0|cjcj|0) = 
1/2 + O (Q), independent of a to lowest order. [\Q\ <C 1 
characterizes the small localized inliomogeneity induced 
by the chemical potential in Eq. (|2.4p .] To capture the 
effects of the lattice, we must cut off the divergence at 
zero argument in the continuum bosonization approxi- 
mation to the lattice correlation function. We do this by 
incorporating a finite range C associated with the nearest- 
neighbor density-density interactions in H^''^ [Eq. (|2.3ap ]. 



C(x„x,0=(0|ct(x,)c(x')|0) 



ctvQ; 



(xj - x'.y + 



(t/2 



7r(a;j - x]) 

kpixj ~ x'j) + t: I dypo{y) 



X sm 



(3.46) 



where po(2/) denotes the initial density profile [Eq. p.BSp ]. 
cn is the normalization constant from Eq. (j3.20p . and 
kp = 7r/2 is the Fermi wavevcctor at half-filling. In 
Eq. p.35|) . K and u denote the Luttinger parameter and 
the sound velocity [for which we will employ the Bethe 
ansatz results in Eq. (|2.5p ]. To compare to the lattice 
quench, we use Eq. (j2.7p for the exponent 17(7). 

The correlator in Eq. p.46p depends upon two length 
scales a and C not defined in the lattice theory. While the 
pure sine-Gordon model results from the limit ^ — ^ 0, the 
parameter a is always non-zero [c.f. Eq. p.20[) ]: its eval- 
uation in the context of the lattice model would require a 
Bethe ansatz calculation.— In comparing to numerics, we 
will fix C = a = 1 (a denotes the lattice spacing), consis- 
tent with nearest-neighbor interactions, but we will treat 
q; as a fitting parameter. Our choices for C and a will not 
prescribe the value C{xj,Xj) = 1/2 + O {Q), except for 
the non-interacting quench cr = 0; rather, we adjust a to 
fit the long-range part of the correlator to the lattice nu- 
merics, since the regularized continuum approximation is 
still expected to behave the worst at short distances. 

Incorporating the same C-regularization into the com- 
ponent correlators in Eq. p.20p and using the result in 
Eq. (|3.13p . one can analyze the "ultraviolet regularized" 
version of the sine- Gordon quench studied in the last sec- 
tion. The regularized post-quench density is expressed as 
the integral 



p{t,x) = ^[po{x-t') + po{x + t')] 



CnOL 



cmct 
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71" dz pq{z) 



t' [f - y) [e + (<' - y fT'^ 



dy- 



*' it' + y)[e + it' + y)r^^ 



dyi / dy2 



(2/1 - 2/2) [C^ + {yi - y2)r^' 

&'\t',y,)G^^\t',y2) - 
+ G(3)(i',yi)G(3)(i',y2) , (3.47) 
+ 2G(2)(t',yi)G(2)(t',y2). 



where G^^'^'^^(t, y) denote the unrcgularized continuum 
Green's functions in Eq. p.9p : the primed time t' = vpt 
[Eq. dUl)]. 
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For interacting initial conditions cr > and a not too 
small, the characteristic "s" shape of the supersoliton ap- 
pears in the regularized sine-Gordon quench. The growth 
of the supersoliton amplitude is terminated after a cer- 
tain cutoff time cx A/vf{MC,)'^ (discussed in more de- 
tail below). The regularized supersoliton is depicted in 
Fig. ini Interpreting the sine-Gordon quench as the con- 
tinuum limit of the lattice model version, we therefore 
anticipate the existence of at least three different dynam- 
ical regimes: (1) < t < 1/{vfM), transient regime, 

(2) 1/{vfM) < t < , "universal" supersoliton regime, 

(3) t > t(^, post-cutoff, non- universal regime. These are 
sketched in Fig. \T0\ 

To obtain an estimate for the cutoff time t^;, we an- 
alyze the asymptotic behavior of the density response 
in Eq. p.47p . We can extract the first correction to 
Eq. p.23[) in the intermediate time window 1/vfM <C 
t <^ t(^; the result is (c.f. Appendix [X} 
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+ {x ^ — x}. 



(3.48) 



where the function Fo-(-z) was defined by Eq. p.24p . The 
correction grows as t^^~'^'>/-^, but with a sign opposite to 
the supersoliton. At intermediate times, the dominant 
effect is the suppression of the supersoliton growth. The 
expansion in Eq. p.48|) is a conserving approximation, 
because Fa{z) integrates to zero (ct > 0). 

In the limit cr — )• 0, the third term in Eq. p.48|) van- 
ishes, as expected for the non-interacting quench (which 
is independent of a and <^). By contrast, ignoring the di- 
vergent prefactor we see that the second and third terms 
precisely cancel for cr = 1. This result obtains because 
the prediction of the unregularized sine-Gordon theory 
suffers a UV divergence for ct > 1; a perturbative expan- 
sion about the C = limit does not exist there. 

We define fi as the amplitude ratio of the second and 
third terms in Eq. (|3.48p , evaluated on the lightcone {x — 



Transient "Universal" amplification/ Non-universal 
regime supersoliton regime cutoff regime 



1/{vfM) 
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FIG. 10. Dynamical regimes of the lattice quench between 
H^'^ and H^-^^ in Eq. (|2.3p . based on considerations of the 
regularized sine-Gordon model, as discussed in the text. The 
quench occurs at t = 0. The cutoff time t^^ can be defined by 
either Eq. 1^:1^ or ((53T]) . 



t'). At argument z = 0, F^{z) equals 2"''^^/^ / r[{2 - 
cr)/4], a value close but not equal to its peak magnitude. 
For fixed ratio fi, we then define the cutoff time 
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(3.49) 



As an alternative, we compare the integral of the absolute 
values of the final two terms in Eq. (|3.48p . We define 



dz \F,{z) 



(3.50) 



Eq. p.50|) can be evaluated numerically. Then, we set 
f2 equal to the ratio of the integrated absolute values 
associated with the second and third terms in Eq. p.48p . 
For fixed f2, we have 



.(2) 



VFiMcy 



f2 2(5--)/4f7^r (-a) 



-|2/(1-<T) 
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(3.51) 



(1 2) 

Both ' diverge as a 



uivci-gc u —r 0, as cxpccted for the non- 
interacting quench. Eqs. p.49p and p.51|) arc rapidly de- 
creasing functions of a that become ill-defined as cr 1. 
Unfortunately, both definitions are also strongly sensi- 
tive to the value that we assign to arbitrary ratio fi^2- In 
particular, with cr = 0.4 or 0.7, for which we present lat- 
tice quench data in Sec. II VI the values of t^^''^^ change by 
several orders of magnitude as fi.2 is swept from 0.1 to 1. 
The sensitivity reflects the very slow in time (fractional 
power law) accumulation of the flnal term relative to the 
second in Eq. (|3.48p . Thus, while the definition of a cut- 
off time with a natural scale A/vf{MC,)'^ is conceptually 
useful, it proves difficult to utilize as a practical tool in 
characterizing finite-size numerics. 

In the ultimate long-time limit t t(^, Eq. p.47p has 
the leading asymptotic behavior (c.f. Appendix 1^ 



p{t,x) 
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(3.52) 



In the regularized continuum quench, the supersoliton 
eventually gives way to a pure translation of the initial 
Gaussian, with a reduced amplitude. This is completely 
different from the single particle evolution resulting from 
"relativistic" confinement < MA ^ 1, discussed in 
Sec. nil Al In that case, the non-dispersive part of the am- 
plitude decays to zero in the long time limit [Eq. p.ll[) ]. 
Eq. p.52p is not a conserving approximation for any cr; 
the missing density is distributed in a long tail neglected 
here. In fact, for C = 1 ^md a > 0.64, the amplitude 
in Eq. p.52p is negative for < cr < 1, which applies 
to Fig. [9] This is the case for the cr = 0.7 and cr = 1.0 
quenches discussed in Sec. IIVI although our lattice nu- 
merics are limited to system sizes much too small to reach 
this regime. 
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3. Band curvature lifetime 

The sine- Gordon theory presented in the previous sec- 
tion accounts only for lattice effects on the initial pre- 
quench state, by way of the (^-regularized correlation 
function in Eq. p.46p . This is one ingredient in the post- 
quench evolution of the lattice density in Eq. (|2.12[) : the 
other is the set of Green's functions G'-^'^'^'^ {t, Xi) ob- 
tained by Fourier-transforming Eq. (|2.1ip . Instead, in 
Eqs. ((07)) ._ ([08)) . and ((3?52|) . we have employed the 
continuum G^^'^'^'^ {t, x) defined by Eq. (|3.9p . which as- 
sumes the Lorentz covariant spectrum in Eq. ()3.26p . 

We find that this regularized correlator + contin- 
uum Green's functions approximation proves adequate 
to model most of the lattice quench numerics presented 
in Sec. lIVI However, to characterize the dynamics in the 
limit of very long times (in a correspondingly large sys- 
tem) , we would need to account for the additional effects 
of band curvature. This is of particular importance for 
the interacting quench, which yields the "regularized" su- 
persoliton in Eq. p.48p {t t^) or its ultimate fate as 
the non-dispersing ghost in Eq. (|3.52p (t ^ t^). These 
disturbances propagate at the "speed of light" vp, which 
is replaced by the maximum band velocity Wmax(-^) in 
the lattice model. [wmax(O) = vp', see Sec. lIV B"2] for more 
details.] As a first correction to the continuum dynam- 
ics, we consider the cubic curvature represented by O3 in 
Eq. dMH). 

The lifetime t^ is defined as the interval post-quench 
during which the cubic curvature can be ignored; a crude 
order-of-magnitude estimate is given by 



t3 



A3 
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(3.53) 



where A is the position space width of the initial density 
inhomogeneity. Eq. p.53p follows from the expansion of 
the band dispersion in Eq. (|2.6p about fcmax such that 
v{k 

max ) — '^max 

fSec. lIVB2l) : 
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where 5k = k — k„ 



5k - -a^5k^ 
3 



a is the lattice constant. For 



a Gaussian packet of width A, the characteristic fre- 
quency associated to the cubic term is u, 
Eq. (p33l) with a = 1. 



/A^, giving 



IV. LATTICE QUENCH RESULTS 

In this section, we present numerical results for the 
XXZ chain quench set up in Sec. [Til A chain with N 
(even) sites and periodic boundary conditions is pre- 
pared in the ground state |0) of Eq. (P3a|) . The 
interaction strength 7 is chosen to reside in the XY 
range — 1 < 7 < 0, so that |0) exhibits gapless power- 
law correlations for the lattice fermions. This state is 
evolved forward in time according to ij'-'', Eq. ()2.3bP . 



For the Gaussian initial state inhomogeneity induced by 
^1°^ in Eq. (|2.4p . we calculate the post-quench dynam- 
ics of the density expectation value p{t,x), Eq. (|2.12p . 
In the generic case of the interacting quench (7 7^ 0), 
the required initial state correlation function C{xj,Xj') 
is computed numerically using the density matrix renor- 
malization group (DMRG) technique. All data shown 
are for a system of = 202 sites. 

We compare the numerical results for the lattice 
quench to the regularized continuum sine-Gordon the- 
ory presented in Sec. IIII D 21 That theory is epitomized 
by the continuum approximation to the initial state lat- 
tice correlation function in Eq. p.46p and the density 
expectation in Eq. p.47p . The regularized sine-Gordon 
model contains two length scale parameters a and C that 
are not defined in the corresponding lattice theory. These 
parameters enter via the initial state correlation function 
C{xj,Xji) in Eq. p.46p . The parameter a determines the 
amplitude of this correlator, while ^ acts as an ultravio- 
let cutoff that renders finite the on-site value of C(a;j, Xj). 
For the non-interacting quench (cr = 0), the continuum 
predictions are independent of C and a. 



A. Non-interacting quench 

The special case 7 = yields a free Fermi gas ground 
state of i/'^*^ Both the initial and final Hamiltonians 
are trivially diagonalized, and we solve for the dynamics 
exactly. For this "non-interacting" quench, the initial 
state correlation function was transcribed in Eq. (|2.14p . 
above. 

We first investigate the quench into the gapless XX 
chain, M = in H^^l Since the low-energy field theory 
description of both the initial and final states is a free 
Fermi gas, we refer to this "EG to FG" quench. 

Thus one prepares a density wavepacket at the origin, 
then simply removes the applied potential and tracks the 
resulting dynamics. Numerical results are depicted in 
Figs. [TT] and [T^ for two different values of A. 

The continuum prediction is a pure translation of half 
the initial density profile to the left and to the right, 
at the "speed of light" vp = 2. The right- moving 
part appears as the first term on the right-hand side of 
Eq. (|3.48p . One can see in Figs . [Til and [T^ that the agree- 
ment between the continuum and the lattice quenches is 
very good, and improves with increasing A. The slight 
dispersion seen for A = 4 in Fig. [TT] can likely be at- 
tributed to the deviation of the band spectrum [Eq. (|2.6p 
with M = 0] from linearity at wavenumbers k ^ 1/A 
away from the Fermi wavevector kp ~ 7r/2. As discussed 
in Sec. IIII D 31 we can associate a lifetime ^3 ^ /vp to 
the presence of the cubic non-linearity in the spectrum. 
Then t3{A = 12) - 860, while tsiA = 4) - 32; the latter 
falls midway in the range of times plotted in Fig. [TT1 

We now turn to non-interacting quenches into a gapped 
final Hamiltonian. Here, we quench from a free Fermi 
gas into a band insulator. The periodic potential in 
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FIG. 11. Time slices of a non-interacting quench with a gap- 
less final Hamiltonian [M = in Eq. (|2.3b|l . a "Fermi gas to 
Fermi gas" quench] at times t — 0, 15, 30, 45, and 60; fainter 
(bolder) traces depict earlier (later) times. Blue solid lines 
are the continuum prediction [the first term of Eq. (|3.48p ]. 
and red dashed lines are the result of exact diagonalization 
of the lattice Hamiltonian. The evolution is symmetric about 
Xj = 0. The relevant quench parameters are Q = 0.10, A = 4. 



Eq. (j2.3bp with M allows for backscattcring umk- 
lapp processes, which open up a bandgap with magni- 
tude 4M. To compare to the continuum theory, we would 
like to reach the scaling limit where all relevant length 
scales in the problem greatly exceed the lattice spacing, 
e.g. A ^ a, 1/M ^ a, while keeping A <C Na. In 
addition, we restrict our quench parameters to the "non- 
rclativistic" transport regime MA > 1, as explained in 
Sec. IIII A[ so as to avoid confusing the putative super- 
soliton (in the interacting quench, below) with relativis- 
tic propagation induced by excessive "squeezing" of the 
initial density disturbance relative to the Compton wave- 
length. Specifically, for all data presented subsequently 
we will fix the product MA = 3/2, and examine four 
wavepacket widths A 4, 6, 12, and 20, yielding the 
respective band-gap parameters AI ~ 3/8, 1/4, 1/8, and 
3/40. 

The application of the staggered potential causes ad- 
jacent site occupancies to 'polarize' opposite to one an- 
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FIG. 12. The same as in Fig. [TT] but with A = 12. 



FIG. 13. Time slices of a non-interacting quench into a 
gapped final Hamiltonian (a non-interacting Fermi gas to 
band insulator quench) at times t = 0, 15, 30, 45, and 60; 
fainter (bolder) traces depict earlier (later) times. Blue solid 
lines are the continuum predictions and red dashed lines are 
the result of exact diagonalization of the lattice Hamiltonian. 
The continuum data results from a numerical integration of 
Eq. (|3.47|l with cr = 0. The evolution is symmetric about 
Xj — 0. The relevant quench parameters are Q — 0.10, A = 4, 
M = 3/8. 



other, but this small-scale density effect is not one in 
which we are interested; the staggered potential is merely 
a tool to induce a gap in the spectrum. We henceforth 
present results for the relative particle density, given by 



dp(t,Xj) = p{t,Xj)Q - p{t,Xj)Q^o, 



(4.1) 



i.e. we subtract the time-dependent density profile orig- 
inating from a spatially homogeneous {Q = 0) initial 
state. 

Figs. [2] (in the Introduction) and [T3] show the result- 
ing post-quench dynamics for two of the four (A, M) 
pairs given above. The dynamics are strongly disper- 
sive, in stark contrast to the ultrarelativistic propaga- 
tion seen in the FG to FG quench. We see that the ini- 
tial Gaussian inhomogencity broadens gradually and does 
so more slowly for larger values of the band gap; these 
non-interacting quench dynamics arc grossly similar to 
the non-rclativistic single particle wavepacket depicted 
in Fig. [SJa). The behavior is generic and we find it to 
occur for a wide range of non-interacting quench param- 
eters satisfying the non-rclativistic condition MA > 1. 

The dispersion arises from a combination of the strong 
band curvature near the (non-interacting) Fermi point 
and the relatively weak occupancy of the conduction 
band induced by the quench, as evidenced by the corre- 
sponding distribution function plots in Fig. [5] with a = 0. 
We emphasize however that the global momentum dis- 
tribution in Fig. |S] does not encode information about 
the inhomogencity; for this purpose one should consult 
the Wigner function, as discussed in Sec. IIII C 21 and 
Appendix [Cj For the continuum theory, we find that 
the "local" velocity distribution for the non-relativistic, 
non-interacting quench exhibits a strong suppression of 
velocities v > vp/MA [Eq. ((C7)) in Appendix [C], due 
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to Pauli-blocking in the initial Fermi gas ground state 
CSec. 1111021) . 

The continuum curves in Figs. [2] and [13] obtain from 
the numerical integration of Eq. ()3.47p . with a = 0. The 
agreement between the lattice and continuum quench dy- 
namics is generally excellent. For the A values consid- 
ered, we observe negligible sublattice staggering in the 
lattice Sp(t, xj); such behavior is a good indicator of the 
near complete separation of the smooth and staggered 
components of the density. This is consistent with the 
retention of only the smooth component of the initial 
inliomogeneity in the regularized sine- Gordon theory of 
Sec. 1111 D2t ^ 



B. Interacting quench 

We turn to the most interesting case of an interacting 
initial state, 7 ^ in Eq. p.3a|) . Because the final state 
is still non-interacting, the dynamics are exactly given by 
Eq. (|2.12p above, but the initial state correlation function 
C{xj,Xji) cannot be obtained via elementary means. To 
achieve this task, we employ the density matrix rcnor- 
malization group (DMRG)S,"-ii due to its ability to treat 
relatively large interacting one-dimensional systems. All 
calculations were performed on a chain of size L ~ 202 
(so that the number of fermions at half-filling, L/2, is 
odd) with periodic boundary conditions (PBCs). In stan- 
dard DMRG, the relative error introduced with PBCs is 
significantly larger than that obtained with open bound- 
ary conditions (OBCs). To achieve the relative error ob- 
tained with m states per block using OBCs, one would 
need states when using PBCs. This results in greatly 
increased computational times which scale as mP with 
PBCs as compared to with OBCs. Efficient methods 
to improve DMRG's ability to handle PBCs are still on- 
going topics of research (see e.g. Ref. [t^ and references 
therein). In spite of the above considerations, we found 
the use of PBCs was necessary to mitigate dynamical 
boundary effects appearing during the quench process. 
In all calculations presented, we kept up to 200 states and 
performed eight sweeps in the DMRG algorithm, yield- 
ing a truncation error (discarded weight) on the order of 
10~^. We tested the combination of DMRG and exact 
time evolution for the non-interacting quench by compar- 
ing to the results of exact diagonalization, presented in 
the previous section. 

Although our DMRG calculations are complicated 
by the use of PBCs and a spatially inhomogencous 
Hamiltonian, one could in principle imagine perform- 
ing DMRG calculations for larger systems. Unfortu- 
nately, although the (non-interacting) dynamics are triv- 
ially written down, they suffer from quite poor poly- 
nomial scaling with system size. Namely, forward and 
backward Fourier transforms [each requiring 0{N) oper- 
ations] for each of the N sites out to times scaling with 
the size of the system yields a dynamics algorithm which 
scales as 0(7V^). Calculating the dynamics for systems 
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FIG. 14. Comparison of the initial state correlation function 
C{xj, 0) defined in Eq. (|2.13|l . as predicted by the regularized 
(a = 0.75, C = 1) LL theory [Eq. (|3.46|) . blue dots connected 
by lines] and as calculated with DMRG (red open circles). 
The inset is a close-up of the same data. The relevant param- 
eters are a = 1.0, Q = 0.0. 

much larger than those considered here is currently pro- 
hibitive. 



1. Initial state correlator 

As a first analysis, we consider the initial state corre- 
lation function C{xj,0) [Eq. (j2.13|) ]. LL theory predicts 
interaction-dependent power-law behavior in correlation 
functions. At large separations, Eq. p.46p yields 

C{x,,0)^\x,\-(^-^'\ (4.2) 

where the exponent a is taken as the Bethe ansatz result, 
Eq. (|2.7p . Fig. [14] compares the correlation function cal- 




FIG. 15. Envelopes of the correlation function C(xj, 0) calcu- 
lated by DMRG for the dynamical exponents a = 1.0 (black 
circles), 0.7 (red squares), and 0.4 (green diamonds) - all for 
Q = 0.0. The dashed lines are a guide to the eye, exhibiting 
slopes -2, -1.7, and -1.4 (bottom to top). 
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culatcd numerically by DMRG to the regularized contin- 
uum prediction for a = 1.0, with (3 = 0. [For the cases of 
non-zero inhomogencity with small \Q\ considered below, 
the Q- and A-depcndcncies of C(xj^ Xji) are very minor.] 
To fix the continuum result in Eq. p.46p . we make the 
physically motivated choice ^ ~ a ~ 1, associated with 
nearest-neighbor density-density interactions on the lat- 
tice, while wc adjust the scale-setting prcfactor a to best 
match the DMRG calculated correlation function at large 
separations. This approach yields cr-dependent values of 
a, specifically a = 0.75, 0.64, and 0.50 for a = 1.0, 0.7, 
and 0.4, respectively. The agreement between the lattice 
and continuum correlation functions is seen in Fig. ll4l to 
be excellent after such a fitting procedure; similar agree- 
ment is obtained for the other values of a. The regu- 
larization parameters obtained in this manner are also 
employed in the subsequent continuum calculation of the 
interacting quench dynamics. 

The power-law prediction emerging from LL theory is 
seen to be very robust. In Fig. [151 '^e plot (in log-log 
scale) the envelope of the DMRG correlation function for 
the interaction strengths yielding exponents a = 1.0, 0.7, 
and 0.4. The deviations occuring at the largest separa- 
tions are an artifact of the numerics. 
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FIG. 16. Time slices of an interacting quench with a — 
1.0 into a gapped final Hamiltonian (an interacting Luttinger 
liquid to band insulator quench) at times t = 0, 15, 30, 45, and 
60; fainter (bolder) traces depict earlier (later) times. Blue 
solid lines are the continuum predictions and red dashed lines 
are the results of DMRG calculations combined with exact 
time evolution of the lattice Hamiltonian. The continuum 
data obtains from a numerical integration of Eq. p.47|l , with 
a = 0.75 and C = 1- The evolution is symmetric about Xj — 
0. The relevant quench parameters are Q = 0.10, A = 4, 
M = 3/8. 



2. Maximum band velocity 

Below we compute the quench dynamics originating 
from an interacting initial state characterized by the cor- 
relation function analyzed above. We first pause to dis- 
cuss a time-rescaling procedure adopted in the follow- 
ing. In Sec. nil C 11 we considered the static post-qucncli 
distribution functions n±{k) for conduction and valence 
band fermions [excitations of the final band insulating 
Hamiltonian H'^^^ or H'^^^] in the lattice and continuum 
quenches. Continuum results for the non-interacting and 
interacting quenches are given by Eqs. p.3ip and p.33p . 
Fig. IS] shows lattice quench results for n+(fc) obtained 
from the DMRG initial state correlation function associ- 
ated with the four values of M considered in this section. 
Each subplot exhibits traces for cr — 0, 0.4, 0.7, and 1.0. 

Fig. [5] indicates that the final state distribution of ex- 
cited particles for an interacting quench (cr > 0) extends 
deep into the conduction band. Although the final par- 
ticle spectrum in Eq. (|2.6p is quadratic at low energies 
near the band center, the distribution induced by an in- 
teracting quench stretches into the linear "relativistic" 
regime of the spectrum and beyond. In the ungapped 
case, the slope vp ^ 2 ioi k just above kp = 7r/2; with 
M > 0, the maximum group velocity of the band struc- 
ture is modified to Wmax(-^^) < 2. It is this velocity with 
which we henceforth rescale time in the continuum cal- 
culations, t' = Wmax(Af)i in Eq. ((XTf)) . 

For Ek in Eq. the velocity v{k; M) = dEk{M)/dk 

evaluates to 

2 sin(fc) cos(fc) 



This equation is maximized at a wavevector fc„iax satis- 
fying 

cos(fc„,ax) = \] + ~m'\ (4.4) 

Inserting the solution of Eq. (|4.4I) into Eq. (|4.3|) yields the 
maximal band velocity Wmax(-^)- In Fig. |8l the position 
of fcniax for each value of M is indicated by a dashed 
vertical line. 

The velocity rescaling procedure outlined above was 
not adopted in the continuum non-interacting quench 
data exhibited in Figs. [2] and [131 For the case tr = 0, 
Fig. m indicates that the linear regime is only weakly pop- 
ulated for all but the largest value oi M = 3/8 considered 
here; see also Eq. p.31|) . As a consequence, the strongly 
dispersive dynamics in the "non-relativistic" transport 
regime MA > 1 are dominated by the low-fc bandstruc- 
ture. This picture is confirmed by the excellent agree- 
ment between lattice and continuum results in Figs. [H 
and[T31 and by the local velocity distribution obtained for 
the non- interacting, non-relativistic quench in Eq. (|C7p . 
For the interacting quenches considered below, the rescal- 
ing of the velocity is not a systematic incorporation of 
bandstructure effects into the continuum Green's func- 
tions defined by Eq. (|3.9p : aspects of ultraviolet band 
curvature beyond the linear regime have been neglected. 
Band curvature effects are expected to become important 
at post-quench times t later than t^, defined as the cubic 
dispersion lifetime via Eq. (|3.53p . 



3. Coherent relativistic wave propagation: 
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FIG. 17. Interacting quench with a = 1.0 as in Fig. 1161 but 
with A = 6, M = 1/4. The continuum result [numerical inte- 
gration of Eq. (|3.47|) ] can be obtained for much larger times 
and system sizes than is currently practical with the interact- 
ing numerics (DMRG-l-dynamics). The top panel shows the 
continuum evolution over a window of length 600. 
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FIG. 18. Interacting quench with a = 1.0 as in Fig. 1161 but 
at times t = 0, 12, 24, 36, and 48, with A = 12, M = 1/8. 



the "regularized" supersoliton 



In Figs. [T6ll24l we present the interacting quench dy- 
namics associated with three different vahies of the ex- 
ponent defined by Eq. (|2.7p . which characterizes 
the initial interacting spin chain described by H^''\ The 
values we choose are a = 1.0 (Figs. [T6HT9)) . cr ~ 0.7 
[Figs. EOl Ell and [1] (in the Introduction)], and a = 0.4 
(Figs. [22H111); these respectively correspond to interac- 
tion strengths 7 = —0.913, —0.872, and —0.790, receding 
from the ferromagnetic transition at 7 = — 1 (Fig.|4]). In 
the unregularized, pure sine-Gordon theory described in 
Ref . I43I and in Sec. lIIIBl the quench yields the prediction 
of the supcrsohton for < cr < 1 [Eq. (|3.23p and Fig. [6] , 
while cr = 1 marks the onset of an ultraviolet divergence 
that must be regularized, as in Eq. p.46p . 

We show data for A = 4, 6, 12, and 20 in Figs. [H- 
[T^ respectively; except for A = 20, the same values 
appear in Figs. [T] and [^OH^ We emphasize that all 
quenches have M = 3/2 A, the same relationship imposed 
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FIG. 19. Interacting quench with a = 1.0 as in Fig. 1161 but 
at times t = 0, 10, 20, 30, and 40, with A = 20, M = 3/40. 
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FIG. 20. Time slices of an interacting quench with a = 0.7 at 
times t = 0, 15, 30, 45, and 60. The continuum data obtains 
from a numerical integration of Eq. (|3.47p . with a = 0.64 and 
= 1. The relevant quench parameters are Q = 0.10, A = 4, 
M = 3/8. 



for the non-interacting case. This constraint puts all of 
our quenches in the "non-relativistic" transport regime, 
as discussed in Sec. IIII Al A single particle wavepackct 
with MA = 3/2 shows only slow broadening, similar to 
the non-interacting quench data in Figs. El and 1131 

The difference in density dynamics for the interact- 
ing quenches shown in Figs. [T] and [T6ll24l as compared to 
the non-interacting versions in Figs. and 1131 is remark- 
able. For all interacting parameter sets investigated, we 
observe a strong separation of dispersive dynamics local- 
ized near the origin (the center of the chain and of the 
initial Gaussian inhomogeneity) , and well-defined left- 
and right-moving wavepackets that propagate away from 
the origin showing minimal dispersion in their spatial ex- 
tents. [Only the right-moving packet is depicted; the left- 
mover is an exact mirror image for the Gaussian initial 
condition in Eq. (|2.4p .] Furthermore, we find empirically 
that these wavepackets travel "relativistically," i.e. at the 
maximal band velocity Vjna.x{M) determined above. By 
tracking the peak of the right- moving wavepacket, we are 
able to extract its propagation speed, which we plot in 
Fig. [25] on top of Wmax(Af). The error bars shown there 
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FIG. 21. Interacting quench with a = 0.7 as in Fig. [20l 
but with A = 6, Af = 1/4. The curve marked "asymptotic" 
is the analytical result for the "regularized supersoliton" in 
Eq. H3.48|l . The top panel shows the numerical continuum 
evolution over a window of length 600. 



FIG. 23. Interacting quench with a — 0.4 as in Fig. 1221 
but with A = 6, M = 1/4. The curve marked "asymptotic" 
is the analytical result for the "regularized supersoliton" in 
Eq. (|3.48|l . The top panel shows the numerical continuum 
evolution over a window of length 600. 
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FIG. 22. Time slices of an interacting quench with a — 0.4 at 
times t = 0, 15, 30, 45, and 60. The continuum data obtains 
from a numerical integration of Eq. (|3.47p . with a = 0.50 and 
C — 1. The relevant quench parameters are Q — 0.10, A = 4, 
M = 3/8. 



FIG. 24. Interacting quench with a = 0.4 as in Fig. [221 but 
at times t = 0, 12, 24, 36, and 48, with A = 12, M = 1/8. 
The curve marked "asymptotic" is the analytical result for 
the "regularized supersoliton" in Eq. (|3.48|l . 



originate solely from the linear fit (of peak position vs. 
time) . The deviation seen at low M (wide A) for weak in- 
teraction strengths likely originates from the inaccuracy 
in determining the exact peak location of such a shal- 
low, wide wavepacket as well as from possible transient 
distortion of the wavepacket's shape over its initial time 
evolution. 

The continuum data in Figs. [T61I24I was obtained by in- 
tegrating the regularized sine- Gordon result in Eq. p.47|) 
numerically, using C = 1 ^-nd the values of a quoted in 
the figure captions. Because the velocity renormaliza- 
tion scheme employed docs not represent a fully system- 
atic incorporation of lattice dispersion details into the 
continuum Green's functions in Eq. p.9p . we anticipate 
poorest agreement between lattice and continuum data in 
the dispersive "tail" dynamics occurring near the origin. 
By contrast, we find very good agreement for the propa- 
gating wavepacket's speed and overall shape, despite the 



crude phenomenological regularization of the continuum 
initial state correlator. 

The largest discrepancies between continuum and lat- 
tice predictions occur for the smallest A = 4, and are 
particularly pronounced in the parameter set {a = 1.0, 
A = 4, AI = 3/8), Fig. [H A significant deviation for 
A = 4 is also observed in the non-interacting, M = 
quench shown in Fig. [TT] Taken together, these results 
suggest that band curvature at the ultraviolet scale k 
fcmax(-M) + 1/A becomes important in the lattice quench 
for this case. This behavior is not accounted for in our 
continuum calculations, which instead assume the mas- 
sive Dirac fermion spectrum Sk = Wmax(-^^)Vfc^ + M'^ 
for the final state Hamiltonian H'-^\ Eq. p.ip (here k is 
measured relative to kp). For A = 4, the cubic curvature 
lifetime in Eq. ^ A^/v^^^{M) = 46, within the 

range of plotted time slices in Fig. [121 For A = 6, our 
estimate for ^3 leaps to 140. 

What can we say about the supersoliton identified in 
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FIG. 25. The maximum band velocity versus the band gap 
parameter, M, as calculated by Eqs. 1)4. 3[l and (|4.4|) [black 
dashed line] and the propagation velocity of the relativistically 
moving wavepacket in the interacting lattice quench for the 
values a — 0.4 [green diamonds], 0.7 [red squares], and 1.0 
[black circles]. 




-100 



7 



FIG. 26. Three dimensional view of the Fermi gas to Fermi 
gas lattice quench (cr = 0, M = 0) with Q = 0.1, A = 6. The 
cyan line demarks the maximal propagation speed, Wmax = 
VF = 2. 



the unregularized, continuum sinc-Gordon quench stud- 
ied in Ref. 113, reviewed in Sec. lIIIBt ' The supersoliton is 
defined as the asymptotic, long time (t S> l/vpM) result 
for the pure sine-Gordon model transcribed in Eq. p.23p . 
exhibited in Fig. |6l The supersoliton propagates ultrarel- 
ativistically at the "speed of light" vp^ has a particular, 
non-dispersing "s" shape and an amplitude that grows 
in time as f^^^. As articulated in Sec. IIII the super- 
soliton arises due to quasiparticle fractionalization. Frac- 
tionalization of the initial, interacting LL state relative 
to the gas of propagating post-quench fermions induces 
a power-law excitation of large momenta in the "local" 
Wigner distribution function n{k;R), as exemplified in 
Eq. p.39p . By contrast, in the non- interacting quench 
momenta jfcj 3> 1/A are exponentially suppressed as a 
consequence of Pauli-blocking (Sec. IIII C 2p . Through 
the massive post-quench dispersion, the Wigner function 
translates into a "local" velocity distribution. In the non- 
interacting quench with AIA ^ 1, only small velocities 
V < vp/MA are significantly excited [Eq. (jC7p in Ap- 
pendix [C]; in the interacting case, a non-integrable di- 
vergence appears at the "speed of light" vp [Eq. (jC8a|) ]. 
irrespective of Af A, signaling the presence of the super- 
soliton. Although rcgularization of the LL correlation 
functions C^ix^x') in Eqs. ((X^ or ((XIS)) at short dis- 
tances ultimately cuts off this divergence [Eq. (|C8b[) ]. the 
distinction between the interacting and non-interacting 
quenches survives, because the power-law behavior in 
C^g{x^ x') at large distances is enough to undermine Pauli- 
blocking, for a > 0. 

In the XXZ lattice quench studied here, we do not ob- 
serve the characteristic "s" shape of the supersoliton in 
any of the lattice data. In particular, the density fiuctua- 
tion 5p{t,Xj) appears strictly positive for the interacting 
quenches, although negative excursions are observed for 
non- interacting quenches. Figs. and [T51 In each Fig.lTTl 
mi and[23](<T = 1.0, 0.7, and 0.4, respectively, ah with 
A = 6), the slim upper panel shows the continuum evolu- 



tion [numerical integration of Eq. p.47|) ] over a window 
of length 600, corresponding to a system size 6 times 
larger than that used for the lattice quench. For cr = 1.0 
and 0.7, the continuum data shows the emergence of a 
negative peak at times and positions much larger than 
could be accessed in the numerical lattice study. The 
relatively good agreement between the lattice and con- 
tinuum results over the 100-site windows in Figs. [17] and 
[^suggests the possibility that the "s" -shape can appear 
in the lattice quench, but only larger system size studies 
can resolve this question;^ 

A more precise way to analyze the interacting lattice 
quench data is to employ Eq. p.48p . This equation de- 
scribes a type of "regularized" supersoliton that appears 
at intermediate timescales: the first two terms are the 
pure sine- Gordon theory result, while the third term is 
the first correction due to a non-zero ultraviolet rcgular- 
ization parameter C^. Eq. p.48p obtains from asymptotic 
analysis of the "exact" regularized sine-Gordon result in 
Eq. ([OT]) . vahd for vpM t < and < cr < 1; see 
also Fig. Uni The cutoff time oc A/vp{M()'^ cx A^ 
(since MA — 3/2 here); alternative definitions of the 
cr-dcpendent proportionality constant are provided in 
Eqs. (jOgj) and (|X5T|) . In Figs, g l EH E H andEl we have 
included time series plots of Eq. p.48p for the quench pa- 
rameters transcribed in the captions. For A = 12, Figs.[T] 
and l241 there is rough agreement between the asymptotic 
result, the lattice quench, and the numerical continuum 
integration for the latest time steps plotted. Eq. p.48p 
fails at earlier times, where transient behavior dominates 
both the lattice and continuum. The asymptotic result 
does not fare as well for A = 6, Figs. [21] and [23l al- 
though the lattice data is well-modeled by the numerical 
integration of Eq. p.47p . Since the cutoff time ~ A"^, 
the failure of Eq. (|3.48p for smaller values of A indicates 
the need to retain higher order terms in this expansion; 
the closest agreement between Eq. p.48p and the lattice 
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FIG. 27. Three dimensional view of lattice quenches with 
Q = 0.1 and A = 6, into a gapped final Hamiltonian with 
M = 1/4 for four different interaction strengths, a = 0, 0.4, 
0.7, and 1.0 (top to bottom). The cyan line demarks the 
maximal propagation speed, iimax(Af = 1/4) ~ 1.56. 



and continuum data occurs for intermediate time steps 
in Figs. [21] and [23l Since we have slaved MA = const., 
both tQ due to the UV regularization of the initial state 
and t3 due to cubic curvature in the post-quench Hamil- 
tonian [Eq. p.53p ] scale as A'^, implying that the effects 
of both types of lattice regularization can be simultane- 
ously decreased by increasing A and the system size L 
for fixed L/A. 

We do not observe the amplification of the initial in- 
homogeneity predicted for the pure sine- Gordon case^ 
in the XXZ quenches. The relativistically propagating 
wavepacket produced by the interacting quench shows 
only increasing diminishment of its amplitude, for all pa- 
rameter sets considered. This is consistent with the anal- 
ysis of the regularized sine-Gordon theory in Sec. IIIID 21 
att = tQ, the prcfactor of the second term in Eq. p.48|) is 
proportional to (a/C)"^, which is less than or equal to one 
for the parameters utilized to model the lattice quench. 

In the Introduction, we exhibited in Fig. |3] a series 
of three-dimensional number density evolution plots for 
the lattice quenches with MA = 3/2 and A = 12, with 
cr e {0,0.4,0.7,1.0}. As a benchmark. Fig. \M depicts 
the cr = 0, M = 0, A = 6 "Fermi gas to Fermi gas" 
quench. (In this case, the quench consists merely of 
turning off the initial Gaussian trapping potential in a 
free Fermi gas.) Fig. [^7] is the same as Fig. [31 but for 
the case A = 6. The weak undulations seen in the cen- 
tral peak of Sp(t, x) for the non-interacting quenches (top 
panels) in Figs. 151 and \T7\ occur at the "Zitterbewegung" 
frequency w = 2Mvf = 4M. These oscillations appear 
in the non-relativistic regime for MA not too large, sim- 
ilar to the single particle wavepacket dynamics discussed 
in Sec. [HIS c.f. Eq. (|XTU1) . 



V. SUMMARY AND CONCLUSION 
A. Summary of results 

In this work, we have performed a systematic study of 
spatiotemporal density dynamics in a ID model of lat- 
tice fermions (equivalent to the spin 1/2 XXZ chain), 
following a quantum quench. The ground state of the 
XXZ chain in the gapless XY phase, parameterized by 
the SfS^^i couphng strength I7I < 1, is time-evolved 
by the non-interacting, band insulator Hamiltonian ob- 
tained by setting 7 = 0, whilst simultaneously turning 
on a sublattice staggered magnetic field (chemical po- 
tential). As a probe of the quench dynamics, we intro- 
duced an additional localized inhomogeneity into the spin 
(fcrmion) density of the initial state, and computed the 
subsequent evolution of the density profile expectation 
value p(t,Xj) under the post-quench dynamics generated 
by a translationally-invariant Hamiltonian. 

All quenches studied in this work feature the special 
property that the dynamics are generated by a simple, 
non-interacting band insulator Hamiltonian, character- 
ized by a bandgap AM. By contrast, the pre-quench 
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initial condition is the ground state of a system of inter- 
acting fermions possessing a low energy Luttinger liquid 
description, except for the special case of a free Fermi 
gas with 7 = 0, referred to as the "non- interacting" 
quench. We used the density matrix renormalization 
group (DMRG) to numerically compute the initial state 
correlation function required to determine the quench 
evolution in the interacting case, studying chains of 202 
sites with periodic boundary conditions. 

We identified a qualitative difference in the density dy- 
namics generated by the initial state inhomogeneity for 
the non-interacting versus interacting (7 7^ 0) quenches. 
For an initial state seeded with a Gaussian density bump 
of width A, in the "non-relativistic" transport regime 
(MA > 1) we found only dispersive broadening for 
the non-interacting quench. By contrast, an interacting 
quench with the same value of MA generates coherently 
propagating left- and right-moving density waves, which 
travel "ultrarelativistically" at the maximum band veloc- 
ity i'max(M) of the post-quench spectrum. 

We showed that the lattice quench data obtained here 
could be well-captured by a regularized continuum sine- 
Gordon model. The continuum theory is an ultraviolet- 
modified version of the pure sine-Gordon quench pre- 
viously studied in Ref. |43| . In that work, an ultrarel- 
ativistically propagating density wave dubbed the "su- 
persoliton" was identified as the leading asymptotic con- 
tribution to the exact result for p{t,x), in the case of 
the interacting quench. The supersoliton exhibits a rigid 
shape and an amplitude that grows in time according to 
where <t > (a = 0) for an initial state possessing 
(lacking) interfermion interactions. In the sine-Gordon 
quench, the supersoliton arises due to the relative quasi- 
particle fractionalization of the initial and final (pre- and 
post-quench) Hamiltonians, quantified by the anomalous 
scaling dimension fT/2 of the post-quench fermions in the 
initial state. In this paper, we showed that fractional- 
ization leads to a divergence at vp in the local (Wigner) 
velocity distribution induced by the density inhomogene- 
ity, for the interacting quench. By contrast, in the non- 
interacting (cr = 0), non-relativistic (MA :s> 1) quench 
we demonstrated that Pauli-blocking limits the excita- 
tion to small velocities v < vp/AIA. 

In the interacting lattice quenches studied here, we did 
not observe amplification of the initial density profile, nor 
did we find the characteristic "s" shape of the supersoli- 
ton. We nevertheless established that the propagating 
density waves produced by an interacting lattice quench 
are well-described by the regularized sine-Gordon the- 
ory. For several of the lattice parameter sets studied, 
we demonstrated that the traveling waves of the corre- 
sponding continuum theory do exhibit the characteristic 
supersoliton shape at length and time scales much larger 
than we can access in the lattice version, owing to com- 
putational limitations. We interpret the waves produced 
by the interacting quench as "elementary excitations" 
of the non-equilibrium state; in the XXZ quench, these 
waves are "regularized" supersolitons. Using an appro- 



priate lattice definition for (7(7), we exhibited the strong 
crossover of the post-quench dynamics as a function of 
the interaction strength. 



B. On field theory methods in quantum quenches 

Beyond the particular dynamical phenomena uncov- 
ered in this paper, our work provides additional support 
to the idea that standard quantum field theory tools can 
be useful in studying strongly out-of-equilibrium physics 
in "realistic" microscopic models. This is non-trivial, be- 
cause field theoretic methods are typically employed in 
condensed matter to capture low-energy, long- wavelength 
equilibrium phenomena such as that observed near a 
quantum critical point. In such established settings, lat- 
tice scale details in the form of irrelevant operators are 
often safely ignored, and field theory tools can be used 
to make robust, sometimes even exact predictions, as a 
consequence of universality. On the other hand, a sudden 
quantum quench in a global parameter of a many-particle 
system typically injects an extensive quantity of energy; 
it is not a priori clear that long-wavelength, continuum 
methods can provide a useful description of the resulting 
dynamics. Indeed, in the context of a quench, irrelevant 
operators encoding lattice-scale details present a serious 
formal difficulty: under far-from-equilibrium conditions, 
rcnormalizability of the low-energy field theory is not nec- 
essarily a barrier against their effects. The problem is 
that long-time dynamics can become sensitive to ultravi- 
olet details, even if the renormalizable field theory (i.e., 
the model obtained by discarding all irrelevant opera- 
tors) gives an ultraviolet finite prediction. The difficulty 
is compounded by the fact that operators irrelevant in 
the infrared become relevant in the ultraviolet, rendering 
perturbative treatments useless for long-time predictions. 

In the lattice quench studied in this paper, irrelevant 
operators suppress the amplification effect seen in the 
pure sine-Gordon model, an integrable field theory in 
dimensions. A systematic improvement of the pure sine- 
Gordon theory in order to describe a particular "parent" 
microscopic model would require a non-perturbative re- 
summation of irrelevant operator effects, a difficult task. 
Nevertheless, we have demonstrated that a phcnomeno- 
logical regularization of the sine-Gordon theory (equiva- 
lent to the "rcsummation" of a particular irrelevant op- 
erator, characterizing the finite range C of the nearest- 
neighbor density interactions on the lattice) gives good 
agreement with the XXZ chain dynamics, at least for the 
system sizes considered here. 



C. Open questions and extensions 

A key question is the survivability of the ultrarela- 
tivistic density packet dynamics for longer times and 
larger system sizes. We identified two lifetimes ti^ and 
^3 that characterize the temporal duration, post-quench. 
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over which lattice effects on the dynamics in the ini- 
tial and final Hamiltonians can be safely ignored. For 
fixed MA > 1 (non-relativistic quench), both cx A'^, 
which implies that the lattice effects can be systemati- 
cally reduced by working with larger system sizes L and 
wavepacket widths A, such that L/A is held constant. 

Finally, the effects of interparticle interactions in the 
dynamical evolution pose a particularly interesting ques- 
tion; would the presence of a non-trivial S-matrix for the 
massive, post-quench spectrum of quasiparticle excita- 
tions tend to encourage or retard the formation and/or 
decay of the "regularized" supersoliton? The answer 
likely hinges upon the presence or absence of integra- 
bility for the post-quench Hamiltonian. In particular, it 
would be interesting to study the density dynamics of an 
XXZ chain quench from the XY phase to the gapped, 
Mott-insulating Ising AFM that occurs for 7 > 1. 

The considerable flexibility afforded to tune con- 
trol parameters in ultracold atom experiments, coupled 
with the excellent decoupling of these systems from 
the environment has brought quench physics in (near) 
integrable models within observational reach.i^^ De- 
spite the powerful methods developed to solve equilib- 
rium properties, so far only limited analytical progress 
on non-equilibrium dynamics in integrable models has 
been made,— 1^"— 1^"— with the exception of systems 
that possess an underlying description in terms of free 
particleSi ^'^'^"''!^^'^^!^'^''^^ Numerical work using the time- 
dependent density matrix renormalization group (t- 
DMRG) by Manmana et. al. in Refs. 113 and Hi on XXZ 
chain quenches between and within the XY and Ising 
phases has revealed the "light-cone effect" predicted by 
Cardy and Calabrese^^, as well as evidence for topologi- 
cal defect formatio n^^'^^ upon quenching into the gapped 
Ising phase. These studies were limited to relatively 
small system sizes (50 sites). A variant of t-DMRG was 
used in Ref. [H to investigate the decay of Neel order 
in XXZ quenches, while a hybrid Bethe ansatz/numerics 
approach was used in Ref. |37| to determine the evolution 
of a ferromagnetic domain wall state. Given the complex- 
ity of the pure analytical approaches, it seems likely that 
a numerical (or hybrid) scheme has the best chance of 
addressing the effects of interactions on the post-quench 
dynamics articulated in this paper. 

We emphasize that even in equilibrium, the effects 
of integrability and irrelevant operators on correlation 
functions at non-zero temperature T > remain sub- 
jects of some controversy ^^i^^" — For massive ID systems, 
Sachdev and Damle^i gave well-reasoned arguments that 
transport at T > should be diffusive. This is also the 
naive expectation for a 1+1-D theory, in the absence of 
other special properties. However, Bethe ansatz results 
on integrable models appear to support the possibility of 
a non-zero Drude weight at non-zero temperature, indica- 
tive of ballistic transport.— One might expect that the 
incorporation of an integrability breaking perturbation 
(such as an irrelevant operator) introduces an additional 
timescale, beyond which the space-time retarded Green's 



function for the appropriate observable (e.g. a spin-spin 
correlation function) would transition from ballistic to 
diffusive behavior. 

Lancaster and Mitra— have investigated a quench deep 
into the Mott insulating breather regime of the sine- 
Gordon modelj^ii^ starting from a LL with an inhomo- 
geneous "domain wall" density profile. In this case, the 
post-quench spectrum can be approximated by massive 
free bosons; because there is no fractionalization, the su- 
persoliton does not occur. Quenches of an inhomoge- 
neous LL with a domain wall density profile into the 
breather regime of sine-Gordon, incorporating interac- 
tions, were further investigated in Ref. [3^ using the 
semiclassical truncated Wigner approximation (TWA)<S 
In this case, the authors uncovered a persistent cur- 
rent, which could signal the preservation of ballistic post- 
quench transport. Previous workS^ has shown that the 
TWA provides a good approximation for quenches into 
the breather-dominated regime studied in Ref. Isol . 

With respect to the phenomena discussed in the 
present paper, the TWA is known to fail^ in the "quan- 
tum" (breatherless) regime of the sine-Gordon modelj^i^ 
where fermionic solitons and antisolitons compose the 
spectrum. The supersoliton has been found at the spe- 
cial Luther-Emery point separating the semiclassical and 
quantum regimes, where there are no breathers and the 
fermions do not interact.— In the Ising phase of the XXZ 
chain, there are also no breathers, and the spectrum 
consists solely of interacting, fermionic spinous. More- 
over, the post-quench dynamics of massive, interacting 
fermions may differ between the continuum sine-Gordon 
and lattice XXZ models. 
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Appendix A: Asymptotic analysis 

In this appendix, we sketch the method used to obtain 
the long-time asymptotic results of Eqs. (|3.23p . (|3.48p . 
and p.52p in the text. All three derive from the ex- 
act expression for the "regularized" sine-Gordon quench, 

Eq. dm . 
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All component integrals in Eq. p.47p feature oscil- 
latory Bessel function kernels; these enter through the 
Green's functions G^'^''^'^^t,y), Eq. (|3^ . Defining 7 = 
Mt' , Eq. p.47p can be expressed as 



p{t,x) 



po{x - t') 



2 21'"+ 
+ {x -x}. 



(Al) 



For simplicity, we consider here only the linear response 
to the initial inhomogeneity pi:){x), for the unregularized 
case with C = 0. Then the integrals /1.2 are 



h = 



dz 



(l-z)i 



2(1-|Z|) 



x — zt' 



dypoiy), (A2) 



dZ 
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dy poiv) 
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Jo yi-iz+^Y] Jo [y^l-{Z-^Y] 

(A3) 

The basic method is to slice up the domain of each 
integral into pieces belonging to one of two varieties: type 
(i) regions throughout which one can employ the large 
argument asymptotic series for the Bessel functions, and 
type (ii) crossover domains where one cannot. For type 
(i) regions, the Bessel functions are replaced by cosines; 
in the absence of a point of stationary phase or some 
other obstruction, these integrals can be systematically 
evaluated by repeated integration-by-parts. Successive 
integrations bring inverse powers of 7 = Mt' from the 
cosine argument, which tend to suppress the contribution 
of the remainder in the long time limit. To ensure the 
convergence of the series obtained for a type (i) region, it 
is necessary to carefully consider the specification of its 
boundary. 

Type (ii) regions, as well as points of stationary phase 
appearing in type (i) domains must be isolated and eval- 
uated by expanding the rest of the integrand in the local 
neighborhood. A useful trick to extract the long-time, 
leading asymptotic contributions is to let each region 
boundary vary with 7 according to a power law. For 
example, the dominant contribution to Ii in Eq. ()A2p 
in the limit 7 — > 00 comes from the narrow type (ii) 
region 1 — Szq < z < 1, where < Szq <C 1. We let 
Szo = ; with ip > Q. Then we perform iterated 
integration-by-parts upon the neighboring type (i) region 
with -I + 5z'o < z < 1 - 7""^ (assuming Sz'^ > 0). To 
ensure that this series converges and produces a sublead- 
ing contribution, one leverages an additional constraint 
(an upper bound) upon the exponent for Eq. (|A2p . 
< V < 2 does the job. Knowing the allowed range 



of 1}} in turn determines the character of the type (ii) 
1 — < z <1 integration. In this way, we isolate and 
evaluate the leading contributions to /i,2 in the long time 
limit, obtaining the asymptotic behavior of Eq. (jAip . 

Execution of the above-described program is straight- 
forward, but tedious; details are omitted here. In the 
remainder of this appendix, we indicate the results by 
identifying the key elements leading to the unregularized 
supersoliton formula, Eq. (|3.23p . 

As discussed above, Ii is dominated by the contribu- 
tion near z — \. Expansion of the rest of the integrand 
gives 
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The I2 integration is dominated by the region with 1 — 
SZq < Z < 1, where < SZq <C 1. For the initial density 
profile, we assume the Gaussian bump in Eqs. p.35|) and 
p.22[) . Making the change of variables Z = 1 — u/27^ 
and Zd = ur/j'^, one finds 



I2 ~ -t'j^^-^poix - t') [l2,a + hb] 



where 

l2,a = 
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The dominant contribution to the kernel K,{u) relevant to 
the evaluation of l2b obtains from the large-u behavior of 
the non-oscillatory term 
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Assuming the Gaussian bump in Eqs. p.35|) and (|3.22p . 
Eq. (jA6bp then evaluates to 



h,b ^ 1 
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In this equation, D^{z) denotes the parabohc cyhnder 
function. Combining Eqs. and (|A10|) yields 



2q(Mag)°^(-.) 



X - t' 



, (All) 



where -F'cr(z) was defined by Eq. p.24p . We obtain 
Eq. ((X^ from Eqs. and (|XTT|) . using Ref. IH. 



these read 
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Number charge conservation is associated with the U{1) 
transformation 

1 ~ . n 



Appendix B: Fractionalization in the sine-Gordon 
model 



In this appendix, we demonstrate that the natural 
quasiparticle degrees of freedom in the continuum Lut- 
tinger liquid Hamiltonian defined by Eqs. p.lSp and 
p.l6|) are fractionalized with respect to the tp fermions 
that appear in the continuum insulator Hamiltonian 
Eqs. dSH) and (IXTO)) . 

We begin by defining canonically rescaled boson vari- 
ables 



so that 



$ = VK<j), e = 9/VK, 
so that Eq. p.l6p can be written as 
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On the second line of this equation, we have refermion- 
ized to obtain an expression in terms of some new, effec- 
tively non-interacting Dirac spinor X'^ The field x car- 
ries scaling dimension 1/2 in the Luttinger liquid with 
Luttinger parameter K, and creates or annihilates the 
"natural" propagating quasiparticle degrees of freedom 
in that phase. The x particles propagate at the sound 
velocity u, rather than the bare Fermi velocity vp- 

Comparing Eqs. (|3.15p and (|B2p . we see that the chem- 
ical potential //^"^ (x) has been rescaled by a factor of 
in the x language. This indicates that the x fermion car- 
ries a fraction \/K of the conserved ip fermion number 
charge. We can see this explicitly by considering the 
bosonic expressions for the components of ip and x; in 
terms of the original boson variables (j) and in Eq. (|3.16l) , 



""0, X ■ 
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Finally, we note that x is non-local when expressed 
in terms of V' (and vice- versa) for any K ^ I, since the 
right-hand side of Eq. (jB4|) must then involve a "string" 
in the argument of the exponential, i.e. an integral of the 
Tp current components { J'^, J^} from minus infinity to the 
argument x; see Eq. p.l7p . 



Appendix C: Wigner functions for the post-quench 
quasiparticles 

In this appendix, we define the Wigner functions for 
the particle ak and hole bk operators of the massive, 
post-quench Hamiltonian H^^'' [Eq. (|3.25p ]. in the (reg- 
ularized) continuum sine-Gordon quench. We then tran- 
scribe results for the local velocity "distributions" in- 
duced by the inhomogeneous po{x) [Eq. p.35p ] for the 
non-interacting and interacting quenches. To simplify 
notation we set the Fermi velocity 



1 



(CI) 



in what follows. 

The particle and hole Wigner distribution functions at 
time t = (immediately post-quench) are defined by 



n+{k;R) = J dxd {0\a^ {R ~ ^) a {R + ^) |0), 

(C2a) 

n-{k;R)= J da-de-*^''(0|6l" {R~ ^)b{R+ ^) |0), 

(C2b) 

where |0) denotes the ground state of H^^\ Eq. p.lSp . 
Both the real space density profile po,± (R) (at time t = 0) 
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and the global distribution function n±{k) can be ex- 
tracted from Eq. (|U2l) : 



poAR) 



dk 
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(fc) = y ^ J dRcxp{-iQR)n+{k;R) 

= /^g(0|at(fc-f)a(fc + f)|0). (C3b) 



For a translationally invariant system, the "point-split" 
integration in Eq. (jC3b[) picks up the delta function con- 
tribution at Q = 0; we are to take e — !• at the end of 
the calculation. 

From Eq. p.27p . a and b are related to the right- (-01) 
and left- movers {1P2) via 

a{k) =p{k)i)i{k) - i/3(-fc)V'2(fc), (C4a) 
h{k) ^m4i-k) - if3{-k)4{-k), (C4b) 



where 



e{k) 



(C5) 



and e(fc) = ^/W+M^. 

We define Sn±{k;R) as the linear response to po{x), 
subtracting the homogeneous (global) distribution. Us- 
ing the correlation functions in Eq. (j3.20p and incorpo- 
rating the ultraviolet regularization C, as in Eqs. p.46p 
and (|3.47p . we obtain 



Sn+{k;R) = -5n_{k-R) 
'dq po{q) 



=CNa 



2tt q 



e"^^B{k;q) 



sgn{k + -j Ga 
- sgn [k - 1^ 



k + 



■A 



(C6) 



where 



B{k- q) =/3 (fc _ I) /3 (fc + f ) + /? (f - fc) /? (-f - fc) . 

The kernel ^^(bl; C) is defined by Eq. (|537)) . Eq. is 
identical to the Wigner distribution for the right-mover 
iIji in Eq. p.36|) . except for the M-dcpendcnt "structure 
factor" B{k\q). 

We consider first the non-interacting quench (cr = 0), 
wherein t/o-(b|; C) — n/l. We assume the Gaussian den- 
sity profile po(a;) in Eqs. p.35p and ()3.22p . As dis- 
cussed below Eq. p.38p . for the non-relativistic (A/ A ^ 



1), non- interacting quench, Pauli-blocking slaves the k- 
dependence of Sn^{k; R) to that of the initial density 
profile po{q — '2k), suppressing the contribution of mo- 
menta |fc| > 1/A. Using the dispersion in Eq. p.26p to 
convert momentum to velocity, we obtain the local ve- 
locity "distribution" at i? = (the center of the density 
bump) , 



Sn+{v;R = 0) 



QM 



cxp 



'{MAY 



2v2(A.fA)s 



(C7) 



This equation applies when v 3> 1/MA, for the non- 
relativistic regime (A/A ^ 1) of the non- interacting 
quench (cr = 0). The exponential strongly suppresses 
velocities v ^ 1/MA. 

For the interacting case, we are interested in a "soft 
quench" fSec. lI AT|) . defined as the regime where 1/A <C 
Af i.e. a non-relativistic initial condition, and an 

effective Compton wavelength much larger than the ul- 
traviolet scale C, which is of order the lattice spacing. As 
in Eq. (j3.39p the position and momentum dependencies 
factorize. Converting to velocity, we obtain 



Sn+{v;R) Ci{a)v'' 
Po{R) 



(1 - i;2)1+'t/2 

1 + 2''-^ 



vMC 

Vl - v"^ 



\ 2 / 



• 1-0- ^ 
, 2 ) 



valid for (1/A/A) < w < 1 - {MCf 12, and 



5n+{v;R) C2(cr)W2-i 



Po(-R) (1 - «2)iW4 



exp 



vMQ 



(C8a) 



(C8b) 



valid for 1 — {MQ'^/2 < u < 1. The prefactors in these 
equations are given by 



ci(ct) 



C2(cr) 



(May (tttT 

2'^+ir(i±^ 
/Afa2\"/2 

A C 



r3/2 



21+.T/2 p (its:) ■ 



Eqs. (|C8a|) and (jC8b[) apply to the interacting quench 
with < cr < 1. For <C 1, i.e. a Compton wave- 
length much larger than the lattice spacing, Eq. (jC8a[) 
exhibits a strong non-integrable singularity approaching 
V = 1. For any C > 0, this divergence is ultimately cut 
off, as in Eq. (|C8b[) . The exponential velocity suppression 
in the latter equation is weaker than that in Eq. (jCTp . 
and originates in the ultraviolet behavior of the regular- 
ized Luttinger liquid correlation function, rather than the 
initial density profile. 
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